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ABSTRACT 


This report describes a study task dealing with the structural dynamics 
of solid rockets. The method used in gaining background material and 
the use of this material in evaluating the structural dynamics of the SRB 
are described. Propellant modes which are self-excited are analyzed in 
detail. The necessity of testing propellant modes is established and the 
requirements for this test are described. 
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SYMBOLS AND UNITS 


a = bore radius of propellant core 

A c = area of normal cross-section of motor case (in^) 

2 

A = area of normal cross-section of propellant (in) 

b - outer radius of propellant core (in) 

c a viscous damping constant (# sec/in) or speed of sound in 

chamber gases (in/sec) 

c c = critical damping constant (# sec/in) 

arbitrary constants 

modulus of elasticity of case (psi) 

modulus of elasticity of propellant (psi) 

storage modulus (real part of complex elastic modulus) psi 
loss modulus (imaginary part of complex elastic modulus) psi 
frequency (Hertz) 
allowable = aUowaWe tensile strength (psi) 

= ultimate tensile strength (psi) 

r *i i — > 

= tan 0 (section 4, 3. 2 only) 


E = 
c 

E = 
P 

E 1 » 
E" = 
f = 


g 

G = 
G> = 
G" = 
h 


K ! 

■j 


gravitational constants or [ K '] 1 
shear modulus or modulus of rigidity (psi) 
real part of complex shear modulus (psi) 
imaginary part of complex shear modulus (psi) 
thickness of motor case (in) 




iv 



area moment of inertia of motor case (in^) 



m = 
M = 
M = 
n = 


area moment of inertia of propellant case (in ) 

zero order Bessel function of the first kind 

first order Bessel function of the first kind 

discontinuity factor or spring constant (section 4. 3. 2 only) 

real part of complex spring constant (section 4. 3. 2 only) 

imaginary part of complex spring constant (section 4, 3, 2 only) 

length (in.) 

tangential wave number for acoustic modes or frequency number 
actual mass 

equivalent mass for constant kinetic energy system 
radial wave number for acoustic modes or mode number 


p = pressure (psi) 

q » longitudinal wave number for acoustic modes 

r = radial coordinate 

R q = bore radius of propellant core (in) 

= inside radius of motor case (in) 

R^ = outside radius of motor case(in) 

R = mean radius of motor case (in) 
c 

s, s*= scale factors 

S. F. - safety factor 

t = time 


v 



U T 


u 


u 


z 

W = 

Y o = 

Y = 
1 

Z a 
<5 _ 


radial displacement 
circumferential displacement 
axial displacement 
■W(r) function r only 

zero order Bessel function of the second kind 
first order Bessel function of the second kind 


axial coord in ate 

w ? 1/2 1/2 

” |(1 + g ) -1 j (section 4, 3. 2 only) 


i4 


characteristic value of tranverse wave (m n) 


3 

density (#/ in ) 


a nan 
T = 

Tre" Y rz' Y ez = shear strains 

& = cubical dilation = 3u r /9r 4- u r/r + (1/r) (9 u e /9 e) +3 u z /9 z 

£ = normal strains in the radial, circumferential and axial 

x v z 

Z = <=/c c 

0 = circumferential coordinate 

v = Poisson l s ratio 

P = mass density « Y/g (#sec^/in^) 

a > a >o - normal stresses 

^ 0 Z 


T re’ T rz’ T az 
oi = 


= shear stresses 


l4 m 


,2 ^ 


circular frequency or Ola. < (1 + £ ) 

2 ( 

circular frequency coefficient for the m the mode 


a/2 
i i/2 + ij • 


(section 4, 3. 2 only} 


Laplacian = 9 2 /a r 2 + (l/r) (9 /g r) -f (l/r 2 ) (a 2 /a © 2 ) +3 2 /a z 2 



1.0 INTRODUCTION 


This report describes a technical study- task which deals with the 
structural dynamics of solid rockets. The scope o_f the ta sk w as to 
first determine the dynamic properties of solid propellant such as 
PBAN and then to determine the dynamic ‘ cfiaracterisHpF'o’f a solid 
booster for which analysis and/or test would be required in a solid 
booster stage development program. The' task was performed for the 
Structural Dynamics Section,~Analytical Mechanics Division, Astronautics 
Laboratory, Science and Engineering Directorate, Marshall Space 
"Flight Center, National Aeronautics and Space Administration. 

During the initial phases of this study, a literature search was 
conducted to obtain reference material related to the structural 
dynamics of solid rockets including material properties of the 
propellant. Many applicable papers and handbooks were located 
and are referenced m the bibliography. Those references in the 
bibliography which are of particular interest are abstracted. 

Concurrently with the literature search, inquiries were made to 
U. S. Army and contractor experts to obtain reference material on 
structural dynamics of solid rockets and to determine what, if any, 
type of dynamic problems have been encountered in solid rocket 
development. 

At Redstone Arsenal, U. S. Army Missile Command, the Propulsion 
Directorate (AMSM1-RK) was contacted. It was ascertained that 
there are no case histories of dynamic problems related to Army 
solid rocket development wherein the dynamics during burn are 
involved. Where the rockets are subjected to transportation environ- 
ments for considerable time before use, fatigue m the propellant grain 
design must be considered. 

Two contractors, namely Thiokol Chemical Corporation (TCC) and 
United Technology Center (UTC) were contacted. At TCC, information 
related to experience with Minuteman First Stage and other solid rocket 
motor development was obtained. This information included the following 

(1) Minuteman First Stage test reports 

(2) Several papers related to modal frequency analysis 

(3} Dynamic properties of PBAN (Polybutadiene-acrylic acid- 
acrylonitrile terpolyxner) 
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At UTC information related to the following was obtained: 


(1) Dynamics related to thrust termination 

(2) Effects of propellant on bending modes' 

(3) Dynamic properties of PBAN propellant 

During the study a review of the dynamic propert ies of a sol id 
pr op e 1 1a nt wa s mad eT ~~ 5 mcV SQlfd~'prO'p~ella J nf~ i s a viscoelastic 
material, the theory for a viscoelastic material and its dynamic 
properties are developed in Appendix A, B ecause PR A M i s most 
apt to be used in futu re development of the solid motor for the 
Space Shuttle, i ts p r op er f i"e s" are s tudl ed and are presented in 
Section Z. 3. ~ 

A review was made of various modes of vibration that can exist 
for a solid propellant rocket and are discussed in Section 3, Of 
all the various vibration modes studied, the study of the Fundamental 
Axisymmetric Longitudinal Thickness Shear Modes received the most 
attention. The reasons for interest in this mode are as follows: 

(1) It may occur at approximately 18 Hertz’ which is in the 
frequency range of interest (0-30 Hz). 

(2) It depends only on the density, dynamic shear modulus and 
the inside-outside diameters of the propellant core, 

(3) It can be excited by the longitudinal (closed-closed) motor 
chamber acoustic mode. 

(4) The dynamic shear modulus varies with temperature, pressure, 
frequency and strain level. Specific shear modulus data for the 
service conditions of the rocket motor to be analyzed is therefore 
necessary. 

Extensive analysis of the longitudinal thickness sh ear .m od-e— i-s-d-eis.cribed^ 
in Section 4. Both theoretical solutions and-'NTfSTRA N computer model 
solutions were obtained to show effects of hd-glT'a-nd KwTmadulxrs7~fiig£rand 
strain and real vs complex modulus. Analysis was also performed to 
evaluate length and size effects on a test segment. 

Dynamic testing of the total Shuttle vehicle in the lift-off configuration 
is considered a necessary requirement to evaluate the complex total 
vehicle dynamic analysis. Dynamic testing of the individual SRB stage 
for those modes involving the SRB structure would therefore not be 
considered practical because of the difficulty of simulating the motion 
restraints which are introduced by the External Tank and Orbiter. 
Propellant modes, namely the Lobar and Longitudinal Thickness Shear 
Modes, would not however be affected by these restraints and could be 



evaluated by individual SRB Staee test* 

modes could be evaluate/! K a-o D * Furt kermore,,both of these 
Of the Fundamentir^L" l a t SRB x Segment test ' Because excitation 
by the motor chamber closed'clos ° nEltudlnal Thickness Shear Mode 

because of the difficulty of ilntifylnAhTd™ u SXpeCtei ’ and als ° 

discussed in (4) above seement Z 1 v y mic shear modulus as 

Asymmetric L„gtodi„a!Tuck„I s Shear® M d 

necessary. R e q uirements for th C. 

fr ° m excitatio ” 

of vibration, its frequencv and th L ° ngltudmal Thickness Shear Mode 

to predict precisely The feasibiHtt- eT ^ ° CCUrrence are difficult 
ment for a test to verifv anal • F a Se S ment test and the require- 

It is therefore recommended Sat I segment^' £ fd StUdy ' 
as possible in the SRB development program ^"ducted as early 
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2. 0 PHYSICAL DESCRIPTION OF SRB 


2. 1 Introduction 


The particular solid rocket motor studied is the Solid Rocket Booster 
(SRB) proposed for the Space Shuttle. The physical size and other 
parameters for this booster were obtained from MSFC-S&E-ASTN-PE. 
At the time the SRB parameters were obtained, the development was 
in the conceptual design phase and therefore subject to change. The 
design studied consisted of a high strength steel case with an outside 
diameter of 142 inches and a length of 1325 inches. The rocket will be 
composed from three to five segments. The propellant, with a circular 
bore of 40 inch diameter, will be poured into each segment. The 
internal pressure will be 1000 psi during combustion. 

2.2 SRB Cast Definition 


The case material is high strength steel, D6AC, which has an ultimate 
strength of 195, 000 psi. Using a safety factor of 1. 4 and a discontinuity 
factor of 1. 2 the wall thickness can be calculated using the following 
equations: 

allowable = Ffu (1) 

sTf. 

and Ft allowable = K p R x (2) 

h 

Combining equations the wall thickness is 

h = KpRiS. F. 

F tu 

or h = 0. 6 12 inches 

2. 3 Core Material Properties 

The propellant used m the study is PBAN (Poly butadiene-acrylic acid- 
acrylonitrile terpolymer) which is a viscoelastic material. It contains 
86% solids and 16% aluminum, has an operating temperature range of 
30 to 90°F and density of 0. 064 pounds per cubic inch. The properties 
of viscoelastic materials m general are discussed in Appendix A and 
the properties of PBAN in particular are discussed in this section. 


( 3 ) 

( 4 ) 
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2.3.1 Poisson's Ratio 


In Section 2. 0 of Appendix A, it is shown that the Poisson's ratio 
for viscoelastic materials is 0. 5. This leads to problems when 
using existing computer codes due to the presence of a 1-2 v term 
in the denominator. The handling of this problem varies with 
contractors. Thiokol Chemical Corporation uses a value of 0.495 
and United Technology Center rewrites their computer codes’ to 
accept 0. 5. It is interesting to note that m the development of the 
equations of the thickness shear vibration* mode (Section 4.2) the 
1-2 v term in the denominator drops out. 

2. 3. 2 Plastic Modulus Shear Modulus 

In Section 3. 0 of Appendix A, it is shown that the elastic modulus and 
shear modulus for viscoelastic materials are complex quantities that 
vary with strain, frequency, pressure, and temperature. The effects 
of temperature can be disregarded due to the low thermal conductivity 
and the high burning rate of the propellant. The increase in pressure 
causes an increase in modulus, therefore, since the interest is in low 
frequencies, it will be conservative to use the modulus for atmospheric 
pressure. The effects of strain and frequency should be included in a 
modal analysis of a viscoelastic material. It was shown m Section 3: 0 
of Appendix A that the moduli increase with a decrease in strain. 
During the contacts with the propellant manufacturers the shear modulus 
for both high and low strains were obtained. The shear moduli for low 
strains was obtained from UTC and is shown in Figure 2-1. The 
equation for the shear storage modulus for low strains can be expressed 
as a function of frequency. It is 

G 1 = 2800 f °- 220 _ (5) 

The shear moduli for high strains was obtained from TCC and is shown 
in Figure 2-2. The equation for the shear storage modulus for high 
strain is 


G' = 900 f °' 145 (6) 

It is seen that there is a large spread in the shear moduli for high and 
low strains. It is therefore evident that it is important to determine the 
strain levels in the propellant and the shear modulus at that strain. 
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2.3.3 Damping Ratio 


In Section 3. 0 of Appendix A, it is shown that the damping ratio is 
given by 


? _ _c_ _ _E_^ 

c c 2E' 

The damping ratio for low strain can be found from Figure 2-1. It 
is 0. 18 in the frequency range of interest* (10-50 hertz). The damping 
ratio for high strain can be found from Figure 2-2, It is 0, 12 in the 
frequency range of interest (10-50 hertz). The damping ratio for both 
high and low strains is high, therefore, the propellant modes will be 
hard to excite and the propellant will act as an energy sink. 
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STRUCTURAL DYNAMICS OF A SOLID ROCKET 


3. 0 


3. 1 Introduction 


A general objective for this work was to define the dynamics which are 
related to the interaction between structure and propellant. Since the 
modulus of solid propellant varies with so many parameters, the modes 
which result from resonance within the propellant were considered to be 
those to which this study should focus. These modes are the lobar 
(breathing) and longitudinal thickness shear modes. When considering 
the possibility of coupling of these modes with total Shuttle Vehicle vibration, 
the longitudinal thickness shear is the most significant. Particular- 
emphasis was therefore placed on this mode. During the study, however, 
reference material was obtained on the vibration analysis of the solid 
rocket in general. In order to present this information for future reference, 
other possible modes of vibration will also be briefly discussed in this 
section. 

3. 2 Bending Modes 

Most applications of a solid rocket would require analysis for the 
clamped-free and free-free bending modes. The mode shapes for 
these modes are illustrated m Figure 3-1 for simple beams and can 
be found in most texts on vibrations of beams. The clamped-free 
bending modes would exist prior to release and the free-free bending 
modes could exist after release. The free-free modes of course 
continually change with the decrease m mass as the propellant is 
consumed. 


The pre-release bending modes of the Solid Rocket Booster for the Space 
Shuttle are more complex than the case described above. The two Solid 
Rocket Boosters (SRB) in one concept are attached to the External Tank 
fore and aft at stations 404 and 1515. The forward attachment points can 
transfer direct loads only to the External Tank. The aft attachment points 
are a sliding joint in the longitudinal direction. Direct loads in the radial 
and lateral directions and local bending moments are transferred by these 
joints to the External Tank. Figures 3-2 illustrates a simple model of the 
Space Shuttle for the lift-off configuration. The inflight model would be 
the same except the restraints at the base of the SRBs would be removed. 

No analysis of the Space Shuttle configuration was performed because it 
was not in scope. It can be seen from the above discussion that dynamic 
testing to evaluate the bending modes of the SRB would have to be performed 
on a eotal Solid Rocket Boos ter /External Tank/Orbiter Configuration. 
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Figure 3-2. SRB Lift-Off Model 






It is significant to note that the properties of the propellant have 
little effect on the stiffness of a solid rocket m bending. This can 
be seen by examining the (El) terms for the case and propellant in 
Table 3-1 for the SRB. 


Note that the El term for the total stage cross-section is only 1/2% 
greater than the case itself. Therefore the case stiffness is seen 
to dominate the total rocket stiffness for any analysis. 


Section 

O.D. 

I. D, 

I . 

E 

El 

Case 

142 in 

141 in 

. 5563xl0^in 4 

28xl0^psi 

1. 558x1 0 13 #in 2 

Propellant 

141 m 

40 in 

. 1928xl0 8 in 4 

4000 psi 

. 008x1 0 3 #m 2 

Case St 
Propellant 

— 

— 



1. 566xl0 2 #in 2 


Table 3-1. SRB Stiffness Terms 


The free-free fundamental bending mode for a solid rocket can be 
calculated by a method described m Reference 18. The solution assumes 
a rocket to be a moderately thick beam and therefore includes the effects 
of transverse shear deformation and rotary inertia. The reference 
solves for the fundamental frequency of a solid rocket having similar 
parameters to the SRB. Solution of these equations are very complex 
and therefore the fundamental bending mode of the SRB was not calculated. 
The solution of the example problem in Reference 18 is shown below since 
it illustrates the method used. Figure 3-3 shows the effective geometric 
parameters required in the solution. 


2-2 Cosh yLCos <SL + 
v/here 


Y (x+Y 2 ) 

6 Cx~<S 1 ) 


6 (x-s 2 ) 

y"(x+y 2 ) 


x = r r 2 =e . a , 5 2 =s+a 
c c 


Sinh yL Sin SL = 0 


a = 


oj ■ 
2 


e - £ 


2 (.YAJ T 

<* G c h c 


4 ( YA) rp(J 2 
? E C I C 


4 * 


+ 


(YI 

9* 


fI) T ] 


.X IXyaKj, 

qG c A c 


( Y i ) rp 
gE c I c 


1/2 
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Figure 3-3. Geometric Solid Rocket Parameters 
Effecting the Free-Free Bending Mode 
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Using the following values gives 


R c = 77. 74 in. 

■R l = 77. 474 in. 

R 0 = 24. 73 in. 

h = 0. 526 in. 

A = 2irR c h = 257 in 2 
c 

A p = u(Ri 2 - R 0 2 ) = l6 > 940 in 2 
I c = nR c 3 h = 7. 76 x 10 5 in 4 
I = ir/4 (R X 4 - R 0 4 ) = 2. 8 x 10 7 in 4 
vc = .28#/in 3 
Yp = . 064#/in 3 

( YI)t = Y c Ic + Yp 7 p = 2. 01 x 10^#-m 
( YA) t = Y C A C + Tp Ap = 1. 157 x 10 3 lb/in 
G c = E c /3 = 9.33 x 10 6 in 4 
f = 5. 5 Hz 

where f = u /2 

3.3 Longitudinal Free-Free Rod Modes 

These modes are the longitudinal resonance of the case wherein the 
mass of the propellant must be considered, but the stiffness of the 
propellant is negligible compared to the stiffness of the case. These 
modes would exist only during transportation and/or handling operations 
of the SRB because in the assembled Space Shuttle configuration the 
case is restrained at its attachment points to the External Tank. This 
mode was determined from the following equation as obtained from 
Reference 18. 

£ - n j A c Eg 

2L A c y c + Ap Yp 
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Using the following values representative of the SRB give 
n = 1 for first mode 
A c = 222. 3 in 2 
Ap = 14357. 9 in 2 
P c = .28#/in3 
P p = . 064#/in 3 
E = 28xl0 6 #/in 2 
g = 386 in/sec 2 
L = 1200 in 


fl = 20. 6 Hz 
3.4 

This mode is characterized by a cylinder whose radius is varying as 
some function of time. Reference 18 indicates that the frequency of 
this mode can be obtained from the equation of motion for a differential 
element, d 0, of a solid motor cross-section which is 

dmu + ku = £ (p) 


Case Radial Hoop Mode (Extensional) 


The solution to this equation is given in Reference 18 as follows: 
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Using the following values representative of the SRB gives 


E c 

= 

28 x 10^ psi 

h 

= 

0. 612 in 

g 

- 

386 in/sec^ 

v c 

= 

. 3 

R o 

= 

20 in 

R 1 

- 

70, 388 in 

R c 

= 

70. 69 in 

Yp 

= 

. 064# /in 3 

Yc 

- 

. 28 #/in 3 

m/m = 

1800 


f s 63. 6 Hz 

3. 5 Lobar Modes (Inextensional) 

These inodes are characterized by displacements developed in the propellant 
when the motor case undergoes arbitrary inextensional circumferential 
, deformations. Strain variations within the propellant occur only within the 
lateral plane. The mode shapes of the case and propellant for the first 3 
modes are as shown in Figure 3-2. Methods for calculating the lobar fre- 
quencies are described in Reference 14. From "Frequency vs. Internal 
Diameter to External Diameter Ratios" plots in reference document the 
fundamental lobar mode for the SRB is expected to be approximately 7 Hertz. 

3. 6 Thickness Shear Modes 

These modes are characterized by displacements developed in the propellant 
in the longitudinal direction. The mode shapes for the first 3 modes for 
the axisymmetric and antisymmetric modes are shown for a slice of a solid 
rocket motor in Figure 3-3. The resonant frequency for the longitudinal 
thickness shear modes are circulated below by equations derived by 
Baltrukonis and Gottenburg (Reference 6) for an infinitely long hollow cylinder 
clamped at the outer diameter and free at the inner diameter. The infinitely 
long cylinder condition is one where only longitudinal displacements occur 

(i. e. , no radial or circumferential displacements). It is shown in Section 
4 that the infinitely long cylinder condition is approached by even segment 

lengths of the SRB. 
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Figure 3-3. Longitudinal Thickness Shear Mode Shapes 
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Axisymmetric modes 

£ m = 1/2 it Vg/ p Ri 2 

Using the following values representative of the SRB give 
G = 900 f°- 145 (from Section 2) 
p = . 064/ g 

an( ^ ^ 3 = 2. 74, 6, 85 and 11. 18 (from Reference 6) 
for R 0 /R 1 = 20/70 = .2857 

fj, £2 and fg = 17. 9, 48. 1 and 81. 5 Hz 

Antisymmetric modes 

f lm = 1/2,,!! m Vg/pS ! 2 

Using the same values for'G, P and R2 as above and the following, gives 

fl. , and m = 3. 48, 7.15 and 11.2 (from Reference 6) 
for R 0 /Rj = 20/70 = . 2857 

ffp f^ an d f|3 = 23. 2, 50. 3 and 81. 6 Hz 
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3. 7 Motor Chamber Acoustic Modes 


Self-generating vibration environments depend on the motor chamber 
acoustic modes. Most of the modes will occur at frequencies too high to 
excite structural vibration modes of a solid rocket except the longitudinal 
chamber fundamental mode. Solid rocket manufactures have computer 
programs to accurately compute acoustic modes for even the most complex 
chamber geometries. Some of the techniques are described in References 
29 and 49. 


The fundamental chamber mode for the SRB can be calculated by treating 
the chamber as a cylinder (Reference 18). The acoustic frequencies can 
then be calculated from the following equation: 


m’ n J 


!(S) 2 + (a m n ) 2 
L R„ 


Using the following values gives the frequencies listed in Table 3-3 
below. 


L = 42, 275 in/ sec 

L = 1190 to 1400 (effective length including effects of radial slots 
in propellant and segment joints) 


R = 20 in 
o 

„ = values were taken from Reference 18 and are listed in 

nmn 

Table 3-2 below. 


Nw fl 
m x. 

0 

1 

2 

0 

0 

1.2197 

2. 2331 

1 

0. 5861 

1. 697 

2. 714 

2 

0. 9722 

2. 1346 

3. 1734 


Table 3-2. Characteristic Value ( a mri ) for Transverse Waves 
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Mode (f m , n ) 

Frequency (Hz) 

Wave Motion 

f 0, 0, 1 

15-18 

Pure longitudinal (closed- closed) 

f i, 0,0 

620 

Pure tangential (fundamental) 

f 2, 0, 0 

1027 

Pure tangential (second harmonic) 

f o, 1, 0 

1289 

Pure radial (fundamental) 

f 0, 2, 0 

2360 

Pure radial (second) 

f i,i,o 

1794 

Coupled transverse 

f 2, 1, 0 

2256 

Coupled transverse 

f l, 2, 0 

2868 

Coupled transverse 

f 

2,2,0 

3354 

Coupled transverse 


Table 3-3, SRB Motor Chamber Acoustic Resonant Frequencies 


The longitudinal (closed-closed) acoustic mode, fg q p consists of an 
oscillating pressure wave within the motor chamber. The magnitude of 
the pressure at the closed ends of motor chamber varies sinusoidally. 

This pressure oscillation will introduce oscillating extensions in the 
motor case which can be expected to excite the longitudinal thickness 
shear mode (using G - 1365 for high strains. High response levels are 
not expected to occur however because (1) the propellant is a heavily 
damped material and (2) the longitudinal shear mode is excited by the 
case which is at a node point. A very high driving point mechanical 
impedance is therefore required to vibrate the case. 

3. 8 Importance of Longitudinal Thickness Shear Mode 

The main effort of this study task was directed towards evaluating 
the fundamental longitudinal thickness shear mode. There are several 
reasons for this. First the fundamental axisymmetric thickness shear 
mode can be excited by the motor chamber longitudinal acoustic mode 
(see paragraphs 3. 6 and 3, 7). A thorough analysis of this mode is therefore 
considered necessary. Second, at the fundamental shear frequency and its 
overtones there exists a strong possibility of coupling with SRB flexural 
and other total vehicle (Shuttle) modes. Third, the requirements for a 

vibration test of the longitudinal thickness shear mode led to further study 
to determine the feasibility of a segment test and an optimization of the size 
of the test segment. 
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LONGITUDINAL THICKNESS SHEAR VIBRATION MODE 


4. 0 


4. I Introduction 


This section will be devoted to the longitudinal thickness shear 
vibration mode of the SRB. First the eigenvalues, eigenvectors, and 
strain distribution for a hollow cylinder, which is traction-free at the inner 
surface, will be determined for real values of the dynamic shear modulus. 

This work will then be used to show the effects on the eigenvalues of the 
dynamic shear modulus for high and low strain values. Next, the eigen- 
value analysis for complex moduli will be compared with eigenvalue analysis 
for real moduli. This analysis will be supported with the comparison of 
the results of the real eigenvalue analysis and complex eigenvalue analysis 
and complex eigenvalue analysis of Nastran. The section will be completed 
with the discussion of the effects of length and diameter of a test segment. 

4. 2 Hollow Circular Bar Clamped Along Its Outer Surface 

An exact solution of the three-dimensional elasticity equations of 
motion in cylindrical coordinates will be obtained for the axi symmetric 
thickness- shear vibration mode of a hollow circular cylinder which is traction 
free at the inner surface and clamped along the outer surface. The solution 
will consist of the eigenvalues, the eigenvectors, and the strain distribution. 
The-solution of this problem was presented by J. H. Baltrukonis and 
W. G. Gottenberg in "Thickness Shear Vibrations of Circular Bars" (Reference 6). 


For the axi symmetric thickness-shear mode, the cylind'er is assumed 
to be infinitely long, since only then is pure thickness- shear vibration possible, 
and displacements occur only in the longitudinal z direction. Therefore, we can 
take the components of displacement to be 


u = u =0 
r © 

(la) 

u = u (r, t) 
z z ' ’ ' 

(lb) 


The three-dimensional elasticity equations of motion, in cylindrical coordinates 
are 
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in which 


and 


(x + 2/i) ¥ “ T 10 + v 17 ’ 


•(2a) 



(2b) 



' (2c) 

. 13 (ru r 1 YBw© ^ 3u* 

r 3r rde ‘ 35' 


(2d) 

, du t _ 3«r 3a, „l/9(r« s ) 

3uA 
33 y 

(2c) 


In view of Equation (1) these equations of motion reduce to the following 
single equation; “ 


! 2 2 

! 9 u , 3u 9 u 

I z , Z _ £_ z 

' 2 T r 9r G . .2 


0r 


Dt 


(3) 


If we now take 


: u = W(r) e 

! Z 


iwi 


(4) 


where 'W(r) is a. function of r only, Equation (3) becomes 


• V + • — + W = 0 ' 

r G 


( 5 ) 


where primes denote differentiation with respect to r. This differential 
equation is the zero-order Bessel's equation whose solution is 


; W<r) = Cj J Q J + C 2 Y o 


(6) 


where Cj and C 2 are arbitrary constants and J Q and Y Q are the zero-order 
Bessel functions of the first and second kinds, respectively. 


The strain- displacement relations in cylindrical coordinates are given 


by; 


r 


9 r 


! 1 j. Jk 

1*8 1 0 8 r 


1 


9u 
z 

8z 


(7a) 

(7b) 

(7c) 
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Pd) 


1 V 


, 3u 

I r 


r0 r D 9 


r 4 9x 


■ 7 rz 


D u 8 u 

r , . 4 

' 

(7e) 

0z Or 


. l 9u z 
0 z r 8 8 


(7f) 


^r8 


In view of Equations (I), Equations (7) reduce to 

C. e r ~ € e = e 2 ~ 7 r 0 = 7 q z - 0 


3 u „ 


'rz d~ 


(8a) 

(8b) 


so that we have only one nonzero component of strain. 

The stress -strain relations in cylindrical coordinates are given by 


“ - * ' 

‘ ““ — * ~ - 


- . - 

- 

4 ; 

”r = 2G 1 

[ 6 r + 

V 

A ) 

* . 

1 - 2v 

r~ " 
\ 

^ = 2G ( 


V 


i 

i 

l-Zv 

i 

» V = 2G 

I * 

(s + 

V 

'*) ' 

1 - Zv 


' T r9’ T rz J t 8z = ^ \ n J 


x9 'rz 10a' 


In view of Equations (8), Equations (9) become 


v = cr n = <r = T 0 =’ r n = ° 
r 8 z r0 8z 


8 u 


rz 


= G 


8 r 


(9a) 

(9b) 

(9c) 

(9d) 

(10a) 

(10b) 


and we see that we have only a single nonzero stress component. By 
Equations (4) and (6), Equation (10b) becomes 


rz 



p* 

f r A 



s 

- -- > 

llpw 2 

C 1 J 1 

2 \ 

P<-> . r V 

( 

r 4 

pw 

G 

L \ 

G j + C 2 Y i 

Cl 

G J 
J \ 


( 11 ) 
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The displacement function. Equation (6), and the strain function. 
Equation 11, were developed for arbitrary boundary conditions. The sub- 
stitution of the appropriate boundary conditions will give the solution of the 
axisymmetric thickness-shear mode of a hollow cylinder, traction free at 
the inner surface and clamped along the outer surface. These boundary 
conditions are 


u 


rz 


■=> 


= 0 


= 0 . 


r=a 


(12a) 


(12b) 


Substitution of Equation (12a) into Equation (4) and Equation (6) gives 
0 = C x J o (b^U a ) + c 2 Y o (b^ P J) (13a) 


Substitution of Equation (12b) into Equation (11) gives 


0 = CjJj 



(13b) 


This system of two homogeneous, linear algebraic equations in the unknown 
constants and can have a solution only if the determinant of the coefficients 
of the unknown vanishes. Therefore, the frequency equation can be written as 


J ° ( w •’ 


(a Jo /) - J, (ajpo ) (1 ^ = 0 



or J o ^m^ Y 1 ^ Q mt, ) “ J i ( Q m “ ) Y 0 ^m^ ~ 


(14a) 

(14b) 


where ^ is a circular frequency coefficient given by 


m 


ft 

m 


X 


"5/ pw m 
= bV~ “ " " 


G 


(15) 


the frequency number jna being appended to identify the zero of the frequency 

Equation (14b) e. g. , ft^ is the smallest value of ft for which Equation (14b) is zero, 
ft 2 I s next higher, etc. The zero of Equation (14b) are tabulated in Table 4-1 
for several inner radius to outer radius (a/b) ratios. 
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• m 

a/b 

1 

2 

3 

0.01 

2.405 

5.522 

8.660 

0.1 

2.44S 

5.726 

9.096 

0.3 

2.786 

6.991 

11.39 

0.333 

2.879 

7.312 

11.93 

0.5 

3.588 

9.604 

15.82 

0.667 

5. 104 

14.28 

23.65 

0.7 

5.619 

15.85- 

26.26 

0.833 

9.776 

28.40 

47.20 

0.9 

16.04 

47.24 

78.61 

0.95 

UBsam 

94.36 

157.1 

0.99 

159.0 

476.1 

793.3 


Table '4-1, Circular Frequency Coefficients 


The eigenvalues can be found, by using Table 4-1 and the following 

(16) 


u 


m 


= Q ) 

m i 


P b 


The eigenvectors can be found by combining Equation (4), Equation (6), 
Equation (13a), and Equation (15), and by letting 


S = 



(17) 


The eigenvectors are given by 


u 

z 


b 


r 

J 

c 

L. 



V s2 ™ 1 

V “J 




(18) 


The strain distribution can be found by combining Equation (11) and Equation (13a). 
The strain distribution is given by 

\ z = -SGS! m [jj (b «m) ^ (B n m)] (19) 

O ft 

Equation (16) and Equation (18) will be used in subsequent paragraphs 
to find the eigenvalues and eigenvectors for different values of shear moduli. 

They will also be used for comparison in the test segment length investigation 
using Nastran. 
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4. 3 Effects of Properties on Eigenvalues 

In this section the effects of eigenvalues will be discussed. The 
discussion will begin with the calculation of eigenvalues for infinite length 
cylinder for high and low strains using the eigenvalue equations derived in 
section 4. 2 and the shear moduli of section 2. 3. 2. Then the effects of 
complex moduli on the eigenvalue will be shown. The disucssion will close 
with the presentation of Nastran results to show the effects of high vs. low 
strains and real vs. complex moduli. 

4. 3. 1 Eigenvalues for High and Low Strains 

The eigenvalues for an infinitely long hollow cylinder traction free 
at the inner surface and clamped at the outer surface were found in section 
2. 3. 2 to be 


W m = ft m G » 

(1) 

where fim is tablulated in Table 4-1 for various a/b ratios. An outer radius 
of 70 inches and an inner radius of 20 inches will be used. Therefore a/b is 
0. 286, and ^ is found to be 2. 74 by interpolating in Table 4-1. In section 
2. 3. 2 the shear storage modulus for high and low strains were given as a 
function of frequency, therefore by substitution into equation (1) we can find 
the real eigenvalues for high and low strains. 

4. 3. 1. 1 Low Strains 

• 

The shear modulus for low strains was found to be 


G = 2800 f°* 22 

(2) 

therefore the eigenvalue for low strain is 

= 2. 74. I 2800 f i°' “ x 386 

0.064 x 702 

(3) 

Since 2 ^f = os , 

(4) 

f l = ^3446 f °* 22 , 

X 2 TT X 1 

(5) 

, ,,,. 0.11 

or f^ = 25. 6 fj. 

(6) 
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Solving equation (6) 

f x = 38. 1 hz . (7) 

This corresponds to a shear modulus of 

G f = 6244 psi (8) 

4. 3. 1.2 High Strains 

The shear modulus for high strains was found to be 

G* * 900 f °* 145 , (9) 


therefore the eigenvalue for high strains is 


or 


wj =* 2. 7 


274 
1 “ 2 ir 


, 0.145 

[900 fi x 386 


0. 064 x 70 2 


1108 f. 


0.145 


Therefore fj = 14. 5 f°* 0725 . 
Solving equation (12) 


fj = 17. 9 hz 

This corresponds to a shear modulus of 
G* = 1365 psi . 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 


4. 3. 1. 3 Normalized Eigenvectors 


The eigenvectors for an infintely long hollow cylinder free at the 
inner surface and clamped at the outer surface are given by equation (18) in 
section 4. 2. This equation will now be solved for a cylinder with an inner 
radius of 20 inches and an outer radius of 70 inches. The circular frequency 
coefficient, fi is therefore 2. 74. The eigenvectors are given by 


u z = 70S 


[j 0 (0. 039 r) - 


J Q (2. 74) 
Y 0 (2. 74) 


Y q (0. 039r) 


(15) 


Equation (15) is normalized and plotted in Figure 4-1. These eigenvectors 
will be used for comparison of the computer results. 
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NORMALIZE DISPLACEMENT {u z /u 






4.3.2 Eigenvalue Analysis 


In this section a comparison will be made between real eigenvalue 
analysis and complex eigenvalue analysis for a multi-degree-of-freedom 
using matrix notation. A complex spring stiffness will be used for damping 
in the complex eigenvalue analysis. 


4. 3. 2. 1 Real Eigenvalue Analysis 


The differential equation for an undamped multi-degree-of-freedom 
system using matrix notation is 

[m]|x| +[k]|x| = 0. (1) 

assuming that 

{x| * {a} e P t , (2) 

Equation (1) can be written as 

? 2 [M] + [K] = 0. (3) 


Therefore the eigenvalue is 



(4) 


and the motion is harmonic with a frequency ofW[Kf] / [m] radians per second. 
4. 3 . 2 . 2 Complex Eigenvalue Analysis 


The differential equation for a structurally damped multi-degree-of-freedom 
system using matrix notation is 


[M] {x 

} + jjK* + i K"J = 0. 

(5) 

Assuming that 



I 3 

={a} 

(6) 

Equation (5) can be written as 


s 2 

[m] + [k 1 + iK n ] = 0. 

(7) 
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The eigenvalue can then be written as 

1/2 

KM i (K>'] 


S = 1 


[M ] + [M] 


( 8 ) 


Letting 


g = 


[M 
ET . 


Equation (8) can be written as 

j,,. ) 1/2 

s = p < 1 + ig V 


(9) 


( 10 ) 


Now to take a root of a complex number we must first express it in exponent 
form. Therefore 


jl +ig[ = ll + g 2 | 1/2 e ie 


where tan © = g. 


1/2 


now 


i + igj = .|i + g 2 | 


(ID 

( 12 ) 

(13) 


The exponential form of a complex number may be expanded as 

(14) 


1 7 • . e 

e = cos o' + i sin — , 
& -2, 


which can be written as 

if / xl/2 

e^ = [ 1 + cos e 


/ V'* / \* 

i ~ i 1 4- cos e \ + if l~cos e \ 

\ 2 j \ 2 ) 


x l/2 


Since 


cose = 


:2\l/ 2 * 


(15) 


( 16 ) 


(1 + g T 

Equation (15) can be reduced to 



+ i If 1 . +..gV- 2 / -lL (17) 

2 1 2 (1+g 2 ) T ^ Z J L ^Z 


1/2 


Combining Equation (10, Equation (13), and Equation (17), the eigenvalues can 
be written as 

/ 

72 


s - 


vr 


[( 1 + g 2 ) l/2 + l] 1/2 +i (l+g 2 ) l/2 -l 


(18) 
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If we take the undamped natural frequency as 




i£‘JM3 * 

Equation (18) can be written as 

3 = - — |{1 + g 2 ) l/2 ~lj 


( 19 ) 


1/2 + i up 


H 


a +g^ 


V^ 2 + il 


Therefore the equations of motion can be written as 

(- a+ i co)t 


lx) = |A) e 


where 


and 


a -5*. {U +gV /2 -l} 

vr 


1/2 


“ {(i +g 2 ) 1/2 + li 


1/2 


( 21 ) 

(22) 

(23) 


1/2 


( 20 ) 


Therefore the motion is harmonic with a frequency of oi and decaying at the 
rate a . It is interesting to note that since g is a positive number 


{(1 hg 2 ) 1/2 + l | 1/2 > 2 


(24) 


and the harmonic frequency of a system with structural damping is greater 
than that of the undamped system. 


4. 3. 3 Effect of Moduli on Nastran Results 


This section will show the effects of the high versus low modulus and real 
versus complex modulus on the eigenvalues and eigenvectors of the Nastran results. 
A comparison of the eigenvalues and eigenvectors will be made with the exact solu- 
tion for an infinitely low elastic cylinder from section 4. 3. 1. The computer model 
used in this comparison is 300 inches in length. A comparison of different 
segment lengths is made in section 4. 4. 1. 

Table 4*>2 shows the effects on the eigenvalue caused by using a real 
and complex shear modulus 



Low Strains 

High Strains 

Theoretical Solution 

38. 1 Hz 

17. 9 Hz 

Real Modulus 

38. 2 

17. 9 

Complex Modulus 

38. 8 

18. 0 


Table 4-2. Effect of Moduli on Eigenvalue 
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associated with high strain and low strains. The shear moduli for the 
real solutions are 6244 psi for low strains and 1365 psi for high strains. 

The shear moduli for the complex solutions are 6244 (1+iO. 384} for low 
strains and 1365 (1+ i 0. 245) for high strains. Table 4-2 shows that there 
is good agreement of the eigenvalues for the theoretical solution and the 
normal eigenvalue analysis. The eigenvectors are compared in section 
4. 4. 1. Table 4t-2 also shows that the eigenvalues for the complex analysis are 
higher than those from the normal analysis. This is expected for the reasons 
discussed in section 4. 3. 2. The increase is greater for the low strains 
because the complex part is a higher percentage of the real part. 

4.4 Effects of Size on the Thickness Shear Mode 


This section will discuss the effects of the segment size on the 
eigenvalues, the eigenvectors, and the strain levels. The effect of segment 
length will be shown by using Nastran models of different lengths. The 
effects of varying the segment inner diameter and outer diameter will be 
discussed by studying the equations developed in section 4. 2. To aid in the 
discussion, the eigenvalue problem of a solid circular cylinder clamped at 
the outer surface will be developed. 

4. 4. 1 Effect of Length 

Three different NASTRAN models were used to study the effects of 
segment length on the thickness shear mode. The models were identical in 
all respects except that one dimension of the elements was adjusted so that 
the model could have a total length of 100 inches, 150 inches, and 300 inches. 
The model is discussed in Appendix B. The model was run at both high and 
low values for the shear modulus using normal mode analysis. Table 4-3 
shows the effects of varying the length on the eigenvalue. Also shown are 
the theoretical .eigenvalues for an infinitely long cylinder which were calculated 
in section 4. 3. 1. The eigenvectors for the various length are shown in Figure 
4-2. The eigenvectors for the infinite length are obtained from section 4. 3. 1. 


LENGTH 

HIGH MODULUS 

LOW MODULUS 

100 Inches 

37. 0 Hz 

17. 3 Hz 

150 

37. 7 

17. 6 

300 

38. 2 

17. 9 


38. 1 

17. 9 


Table 4-3, Eigenvalues Versus Length 
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NORMALIZED DEFLECTION 



70 60 50 40 30 20 

RADIUS (INCHES) 

FIGURE 4-2. REAL EIGENVECTORS OF FIRST THICKNESS SHEAR MODE FOR HIGH SHEAR MODULUS 


The eigenvectors shown are normalized, and, thereto re, are not influenced by 
a change in modulus. 

Table 4 -3 and Figure 4 -2 show that the solution of the eigenvalue 
problem depends on the length of the segment. The solution of the problem 
approaches that of an infinite length as the length is increased. The length of 
the segment was not made greater than 300 inches due to the natural boundary 
conditions made during manufacturing. It assumed that a vehicle segment 
will be approximately 300 inches in length. 

4. 4. 2 Solid Circular Cylinder Clamped at th'e Outer Surface 

The eigenvalue problem for an infinitely long elastic solid circular 
cylinder clamped at the outer surface will be solved to aid in the discussion 
of varying the outer diameter of the test specimen for a scaled-down test. 

The solution of the problem will be identical to that in section 4. 2, only the 
boundary condition will change. As in section 4. 2, the equations of motion can 
be reduced and solved as 

W(r) = Cl J 0 (r^27 + C 2 Y 0 (r-^ET )' (1) 

The boundary condition for this case is 

u 2 | = 0 (2) 

r = b 

In order to avoid infinite displacements at the axis of the bar, must equal 
zero and therefore 

Jo ( n m) = O (3) 

The zero values of pm are given in Table 4-4. 


m 

: 

Q m 

1 

. . , . . .. 

1 2. 405 

2 

! 5.520 

3 

8. 654 

4 

11.79 

5 i 

■ 1 

14. 93 


Table 4-4. Natural Circular Frequency Coefficients for 
a Solid Circular Cylinder Clamped at Its 
Outer Surface 
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The deflection is given by 

Pz a S J G {^m L ) (4) 

b b 

and the strain is given by 

I = _sg j (^j£) (5) 

b 

4. 4. 3 Effects of Size on Test Specimen 

The purpose of testing a segment is to reduce the cost of testing 
the total vehicle or stage, but the test will only be valid as long as the parameters 
that are measured in the test are not affected by the size of the test specimen. 

In section 4. 4. 1 it is shown that to produce the pure thickness shear mode 
the test specimen must be approximately 300 inches long. The effects of 
scaling the diameters is a more complex problem. 

Many things must be considered when determining the effects of 
scaling the diameters. The moduli are functions of frequency. The moduli 
are difficult to measure and, therefore, the frequency between the test specimen 
and the vehicle should not change. Since 

~ ( 1 ) 

N Pb 2 


and m is a function of a/b, as b is decreased a would also have to be decreased 
to keep the eigenvalue constant. The limit would be reached when a reaches zero 
and the cylinder becomes solid. The limit for b can be found by using section4. 4 .2. 
Using the values for low strains and assuming that G is only a function of frequency 


38. 1 x 2 it 


2. 405 


4 


6244 x 38 6 
. 064 x b ^ 


( 2 ) 


and solving, b - 62 inches. Therefore, the smallest diameter test specimen with 
no change in the thickness shear mode would be 124 inches if G was only function 
of frequency. But it is shown in Appendix A and in section 2 that the moduli are 
functions of frequency and strain. 

In section 4. 2, it was shown that the strain is a function of the outer radius, b, 
and the ratio of the inner radius to outer radius. The amount of strain variation 
can be seen by considering the term Yo (ftm) in equation 19 of section 4. 2. At 
certain diameter ratios this term can go to zero. When this happens the strains 
can become infinite. This is only in a theoretical sense, but still, a large varia- 
tion in moduli and, therefore, frequency would occur. Therefore, only the full 
scale diameters should be used in vibration testing. 
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4.5 


The longitudinal thickness shear mode is a function of the modulus, 
density, and the physical configuration. The moduli are functions of frequency, 
therefore the value used to calculate the thickness shear mode must correspond 
to the eigenvalue. The moduli are also functions of strain, therefore the value 
used should correspond to the strains. 

Using the complex modulus causes only a slight increase in the eigen- 
value. This was shown both theoretically and with NASTRAN. The increase 
is slight enough to disregard the complex stiffness and use only a normal mode 
analysis unless very accurate results are needed. 

An operation SRB segment should be used as the test segment to evaluate^ 
the thickness shear mode. If, however, the segment length is considerably longer 
than 300 inches, considerations should be given to reduce the length to 300 inches 
for the purpose of total test program cost reduction. The diameter of the test 
segment can not be sealed because the varying strain will effect the eigenvalues. 
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5. 0 SRB TEST RECOMMENDATIONS 

5. 1 Type of Testing 

The type of testing discussed in this section is development testing. 

The primary purpose of this type of testing is to evaluate the design and/ 
or verify the analysis. It should be distinguished from qualification testing 
which is primarily conducted on production hardware to certify the manu- 
facturing processes. Development tests are considered those tests which 
are conducted early in the program on prototype models before static 
^firing and flight testing is performed. 

5. 2 Requirements for SRB Dynamic Testing 

The precedent of previous vibration testing of solid rockets such as 
Minuteman and Titan suggest the need for full scale and individual stage 
vibration testing to evaluate the various modes of vibration including 
modes predominantly propellant oriented as well as structural modes. 
Although discussions with TCC and XJTC indicated no design problems 
were revealed by development testing, they did afford a verification of 
analysis. 

The Space Shuttle contains an External Lox/LH^ Tank, two Solid Rocket 
Boosters and the Orbiter. It is much more complex structurally than the 
Minuteman or Titan Vehicles, Experience with Saturn V dynamic testing 
has shown the deficiencies of dynamic analysis where a large complex 
launch vehicle is involved. Although the techniques of structural analysis 
have improved greatly from the development periods of the Minuteman, 

Titan and Saturn V Vehicles, certain assumptions and analytical modeling 
simplifications are still required. It is these assumptions and simplifica- 
tions which must still be verified by test. 

5. 3 Lateral Bending and Longitudinal Structral Modes 

Dynamic development testing of the Space Shuttle Vehicle provides the 
ultimate configuration for evaluating the SRB because the exact motion 
restraints of the SRB attachment hardware can be realized. All structural 
oriented modes such as lateral bending and-longitudinal modes can be 
ideally evaluated in the total Space Shuttle Vehicle configuration. Separate 

Separate dynamic testing of the SRB for evaluating lateral bending and longi- 
tudinal modes with reasonable accuracies would require elaborate and costly 
simulation of the External Tank and Orbiter interfaces. These interfaces 
would include local stiffness of attachment hardware and mass property 
simulation of the External Tank and Orbiter. It is doubtful that the cost 
of implementing a separate SRB test would be attractive if development 
costs of interface hardware were considered. 
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Because of the above arguments, it is recommended that the 

SRB lateral bending and longitudinal structural modes be evaluated by 

a dynamic test with the SRB in the complete Space Shuttle configuration. 

5, 4 Propellant Modes 

Propellant modes as discussed in Section 3 are the Lobar and Longitudinal 
Thickness Shear Modes. 

The requirement to verify the analysis bf propellant modes stems from 
the fact that the shear and elastic moduli are complex quantities which 
vary with strain, frequency, pressure and temperature. Solution for 
the propellant modes can, therefore, vary considerably based on the 
shear moduli data used in the analytical model. 

Propellant modes are not dependent on the external motion restraints of 
the SRB. Separate full scale SRB stage testing could therefore be conducted 
to evaluate these modes. Also since it is considered necessary to evaluate 
the propellant modes with active propellant, there would be ignition hazards 
associated with propellant mode testing. It is therefore desirable to limit 
the risk and test the smallest configuration. 

The feasibility of evaluating propellant modes with (1) a segment of 
the SRB and (2) a reduced scale segment were evaluated during this study. 
The results of the evaluation, which are summarized in section 4.4. 1 
and 4. 4, 3, show that a segment test would be a satisfactory test specimen, 
but a reduced scale segment would be unsatisfactory. Furthermore the 
length of the test segment has been established at approximately 300 inches. 
Based upon the above arguments, it is recommended that the Longitudinal 
Thickness Shear Mode analysis be verified by a SRB segment test. 

5. 5 Motor Chamber Acoustic Modes 

/ 

Motor chamber acoustic modes are a source of self-generated vibration. 
Only the Fundamental Closed-Closed Longitudinal mode is low enough in 
frequency to excite structural modes of the SRB. Analytical techniques 
are considered adequate to accurately determine this mode. No develop- 
ment testing is therefore recommended. 

5 . 6' Material Sample Testing 

Reference 9 has shown that the elastic modulus of elastomers can vary 
drastically as the amount of filler material is changed. The properties 
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or the SRB propellant will therefore most likely have to be completely 
evaluated after formulation with regards to the amount of solids in. the 
propellant is specified. Furthermore, the data available from literature 
ar_d from TCC and UTC is not adequate to show the specific effects of 
temperature, pressure, frequency and strain level on PBAN. 

In. order to have the necessary data to accurately predict propellant 
modes, it is recommended the sample testing be performed early in the 
development program so that "Dynamic Shear Modulus vs. Frequency" 
plots can be obtained in the range of 5 to 50 Hertz at several strain levels 
from the lowest to the highest strains expected. It is recommended that 
tine effects of temperature be evaluated at 30, 70 and 100 °F, and the 
effects of pressure be evaluated at barometric pressure and 1000 psia. 

There are many test methods that can be used to evaluate dynamic shear 
modulus. The methods described in Reference 10 are used in the pro- 
pellant industry and will probably give the best results. 
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APPENDIX A 

DYNAMIC PROPERTIES OF VISCOELASTIC MATERIALS 




SYMBOLS 


E 

G 

v = 

K 

F 

a = 
e = 
n = 

T — 

E r - 

E ' 


E" = 
<3 



H 

T 


d 


G 1 

G" 


A 


w = 
t — 
E* 



Elastic Modulus 

Shear Modulus 

Poisson's Ratio 

Bulk Modulus 

Force 

Stress 

Strain 

Viscosity Coefficient 
Relation Time 
Relaxation Modulus 
Elastic Storage Modulus 
Elastic Loss Modulus 
Loss Tangent 
Lamping Ratio 

Distribution of Relaxation Times 
Temperature (Absolute) 

Density- 

Temperature Shift Factor 

Shear Storage Modulus 

Shear Loss Modulus 

Magnitude of Harmonic Function 

Frequency 

Time 

Complex Elastic Modulus Magnitude 

Dashpot Strain 

Spring Strain 
Degrees Kelvin. 
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APPENDIX A 


Dynamic Properties of Viscoelastic Materials 


1.0 INTRODUCTION 


The physical properties of viscoelastic materials required for dynamic 
analysis will be defined and characterized in this appendix. These 
properties are elastic modulus, shear modulus, Poisson's ratio, and 
the damping ratio. Only three of the four parameters are indepen- 
dent, since three of the parameters are related by 

E = 2 (1 + v )G. (1) 

The significance of these properties will be discussed as they relate 
to viscoelastic materials. 


2.0 POISSON'S RATIO 


Poisson's ratio is related to the bulk modulus and shear modulus by 
(Reference 19) 


v 


3K-2G 

■ erKTzrr 


which can be reduced to 


V - 


3 - 2 ^ 

TTST 

K 


( 2 ) 


(3) 


It is customary to assume that viscoelastic materials are incompressible; 
therefore the bulk modulus is infinite (K = « ) and Poisson's ratio is one- 
half ( v ss 0. 5), Equation (1) can then be reduced to this simple relationship 


E = 3G 


(4) 


3.° MODULI 

Since the relationship between the elastic modulus and the shear modulus 
for viscoelastic materials is given by equation (4), only one needs to be 
determined. The elastic modulus will be determined because of the ease 
in visualization. To aid in the development of the elastic modulus for 
viscoelastic materials, the stress-strain relationship for elastic materials 
and viscous materials will be presented. They will then be expanded to 
show the stressr-strain relationship for viscoelastic materials. 
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3. 1 Elastic Materials 


In an elastic material (metals) operating below the elastic limit, stress 
is directly proportional to strain. .This can be modeled by a linear spring 
as shown in Figure A-l. 

LLW 



F(t) 
o (t) = E e (t) 

Figure A-l, Model of Elastic Materials 


It is important to note that the dependent variable c (t) is directly 
proportional to the independent variable without a phase shift during 
dynamic excitation. The response to step stress <r c , is shown in 
Figure A-2, and the response to a step strain, £ 0 , is shown in Figure 


A-3. 


£ 



I * 

Figure A-2. Elastic Response to Step Stress 


cr= E e 

o 


t 


Figure A-3. Elastic Response to Step Strain 


3. 2 Viscous Materials 

In a linear viscous material (fluids) the stress is directly proportional 
to the rate of strain. This can be modeled by a linear dashpot shown in 
Figure A-4, 


i F(t) 

a (t) = n d £ (t) 
dt 

Figure A-4, Model of Viscous Materials 
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The response of a viscous material to harmonic excitation can be 
found by letting 



s (t) = 

A sin to t 

(5) 

then 

lift) = 

cuA cos oj t 

(6) 


dt 



and 

ft- 

II 

V o)A cos a) t 

(7) 

or 

cr (t) = 

7j O' A sin ( wt + 

L) (8) 


By observing equations (5) and (8), two important facts can be seen: 

1) Stress leads strain by 90°, and 

2) Stress is a function of the frequency of excitation 

The response to a step stress is shown in Figure A-5. The response 
to a step strain is an impulse function and is not shown. 



Figure A-5. Viscous Response to Step Stress 
3. 3 Viscoelastic Materials 

As the name implies, viscoelastic materials (polymers) have some 
characteristics of elastic materials and some characteristics of 
viscous materials. Many mathematical models have been developed 
for viscoelastic materials. All of these models are various combina- 
tions of the spring and dashpot elements. Only the simplest model will 
be discussed. . 

Viscoelastic materials can be modeled by the spring and dashpot element 
in series as shown in Figure A~6. The relationship between stress and 
strain can be found by substituting the following identities into the equations 
in Figure A-6j 
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///////// 



e (t) = c s (t) +£#) 

cCt) = a s (t) = a d(t) 

v/here s indicates spring 
d indicates dashpot 


Figure A-6. Simple Model of Viscoelastic Materials 


S s (t) = 


g (t) 

E 


and 



substituting 


dt 


( 9 ) 

( 10 ) 


s(t) = 



dt 


(ID 


or d e(t) = d cr(t) + q-(t) (12) 

dt E dt V 


By defining the relaxation time, as 


Y = 


n 

If 


(13) 


equation (12) can be reduced to the following form to give the relationship 
between stress and strain: 


a (t) = 


■p d j 

2 ’3T 


d 

+ 1 j 


dt 

t ) 


(t> 


(14) 


It can be seen that equation (14) is far more complex than that for 
elastic behavior. Although this math model is very simple, it can be 
used to obtain an understanding of the response of viscoelastic materials 
to various forms of excitation. 
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3. 3. 1 Viscoelastic Response to Step Functions 


The response of a viscoelastic material to a step stress is commonly 
called creep. The strain resulting from a step stress is 



e(t}~ 

a o + 

foo dt 

(15a) 

' 

e' 

' n 


or 

m = 

a o -}• 

cot 

(15b) 


It can be seen that initially the spring elongates an amount equal to 
a Q /E and then the dashpot elongates as a function of time. The response 
is shown in Figure A-7. 



Figure A-7. Viscoelastic Response to a Step Stress 

The response of a viscoelastic material to a step strain is commonly 
called stress relaxation. Since strain is not a function of time for a 
step strain 

0 (16) 


+ ojt) = 0 (17) 

T 

The solution of this differential equation is 

a (t) = a 0 e-t/r (1 §) 

It can be seen that the initial deformation goes into stretching the spring. 
Then the stretched spring exerts a force on the dashpot which elongates 
and allows the spring to contract reducing the force. The response is 
shown in Figure A-8. 


fojt) = 

dt 

therefore 

d q (t) 
'dt 
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-t/x 


Figure A-8. Viscoelastic Response to a Step Strain 

The curve in Figure A-8 is called a relaxation curve. The relaxation 
modulus is defined as 

E r (t) = a = “ t/ T = E e _t/ T (19) 

e e 

It is important to note that the relaxation modulus is a function of time. 

The relaxation modulus will he discussed again later. 

3.3;2 Viscoelastic Response to Harmonic Excitation 

The response of a viscoelastic material to harmonic excitation can give 
us important insight regarding the parameters required for modal analysis. 

Letting E (t) = A sin cot (20) 

in equation (14) yields 


do (t) 4- q (t) = AjjE co s to t 

dt x 


( 21 ) 


The solution of this differential equation is 


o(t) = a) EA (to sin a) t + — cos ut) . . 

7 . , (42) 

1 + W 1 2 

~2 

T 


Since 


cos to fc - sin ( tjjt + ^ ) 

2 


(23) 


or 


0 (t) = oj 2 EA sin fcjt + to'EA t (sin, 0 't + ■£) (24) 

u2 7T _2 ' 


V 


1 + a 
12 


o (t) - to 2 t2ea sin to t + ootEA (sin^ t + 2 ) (25) 


1 -fw" 2 ■? 


1 + (,J 
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Since e (t) = A sin wt, it can be seen that there are two components in 
the relationship between stress and strain. One stress component is in 
phase with strain and the other stress component is 90° out of phase with 
strain. This relationship can be written as 



a (t) = (E - + iE 1 ') ( e (t) 

(26) 

or 

<J(t) = E* + iE" 

(27) 


e (t) 

I * I) 

Where E represents the in phase component and E represents the out 
of phase component. They are given for our simple model as 

e' = to 2 - t 2 E (28) 

1 to * 2 ' T 2 

and E 1 ’ g MT E (29) 

1 + (0 2 T 2 

3. 3. 3 Storage Modulus and Loss Modulus 

E 1 is known as the elastic storage modulus and E M is known as the 
elastic loss modulus. It is important to note that the moduli are 
functions of frequency. The relationships between the components 
are shown in Figure A~9 and Figure A- 10. It should be noted that the 
strain lags the total stress by the phase angle 6 . Since the relationship 
between stress and strain is a complex number, it may also be written 
as 

E* + E"= E re" 6 (30) 

where = |(E ! ) 2 + (E n ) 2 j 1/2 (31) 

and <5 = tan**^ ( 22 ) 

E 



Figure A-9. Relationship Between Stress Components 
and Strain 
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3. 3. 4 Complex Modulus and Loss Tangent 


The vector sum of the storage modulus E) and the loss modulus E n , is 
defined as the complex modulus E This is shown in ‘Figure A-10. 



Figure A-10. "Vector Representation of Complex Modulus 


The parameter <5 is called the loss tangent and is related to the damping 
ratio by 


c ~ ^ - tan 6 - E" 

cc 2 El (33) 

3. 3. 5 Shear Modulus 

It has therefore been shown that the relationship between stress and 
strain far our simple model of viscoelastic materials is a complex 
number which is a function of frequency. Although we elected to use 
the elastic modulus for our development, we could have just as easily 
used the shear modulus. In fact, the relationship between the shear 
modulus and the elastic modulus given in equation (4) holds true for 
the complex moduli. Therefore the shear storage modulus, G 1 , and 
the shear loss modulus, G , are given by 

e’ = 3G* (34) 

and E n = 3G M (35) 

3. 3. 6 Relaxation Modulus and Dynamic Mo duli Determin ation 

Equations (9) thru (29) were developed for a one spring-dashpot system 
with a relaxation time of t . By comparing the mathematical model 
response in Figure A-8 with the experimental response in Figure A- 11, 
it can be seen that the model is not representative of measured data over 
long time periods. A better representation can be found however, by 
using an infinite number of relaxation times. Therefore, the relaxation 
modulus can be expressed as 



Ef (t) = Eje ^ + E 3 e- t/T 3, _ _ 

(36) 

or 

CO 

Er'(t) = lim E ± e“ fc/<r i 

(:?) 


•“ 1*1 
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( 38 ) 


Equation (37) can be written in integral form as 
E r (t) = f BJ T )<»-*/? dt + E m 

Ts» O 

E 00 is added after integration because after an infinite time the 
relaxation modulus ceases to decay without having reached zero. 

E 00 represents the value of the relaxation modulus at that time. 

Since the distribution of relaxation times is very broad, it is convenient 
to use logarithmic notation. 


Letting H (In x ) := t E(t) 


/ CO 

H(ln x )e d(ln t ) -f E a 


lm - - 


(39) 

(40) 


H(ln t) can be calculated a number of ways. A crude approximation 
is (Reference 17 ) 


H (In T ) = 


- d [Er(t)| 
" d hit 


(41) 


A more accurate approximation is (Reference 17) 

H (to T) =irkt [ Er w] + 

Equations (41) and (42) are evaluated under the conditions that t =x 
In order to get an accurate approximate distribution of relaxation times, 
the experimental data should cover at least eight to fifteen decades of 
time. Returning to Equations (28) and (29) for dynamic moduli, they can 
now be expanded to give 



and 


/ co 

H(ln T ) 

£0 

/ CO 

H(lnx ) 


t 2 d ( lnT ) + 

ds 11 d(lnt) 

1 + x 2 



(43) 

(44) 


Therefore, the complex modulus for dynamic analysis can be found by 
measuring the response to a step strain. 


From the stress relaxation curve, the distribution of relaxation times 
(H (lnx) ) can be found. This distribution can then be used in the appro- 
priate equation to determine the dynamic moduli. Computer codes have 
been written to aid in the mathematics. To obtain an accurate distribution 
of relaxation times, the relaxation curve must be known over eight to 
fifteen decades of time. To reduce the test time required to generate 
the relaxation curve and to alleviate the difficulties of obtaining data for 
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short times due to the inability of test equipment to apply an exact step 
strain, a time -temperature superposition principle may be used. This 
principle is explained in detail in Reference (17} and is next summarized. 


The time -temperature superposition principle states that in viscoelastic 
materials time and temperature are equivalent to the extent that data at 
one temperature can be superimposed upon data taken at a different 
temperature merely by shifting curves. The procedures involve shifting 
of stress -relaxation curves taken at different temperatures along the log 
time axis until portions of the curves all superimpose to form one contin- 
uous curve. Before the experimental relaxation curves can be shifted 
to make the master curve the modulus values should be corrected for 
density and temperature. This correction is 

E^t) reduced - Er (t) (45) 


where the subscript , o’ indicates a reference temperature. After this 
correction has been made, the curves may be plotted and then shifted to 
obtain a master relaxation curve. The amount of shift is given by the 
following equation. 


log A t = log t/t Q 


-8. 86 (T-243) 
lOXTb T-243 


(46) 


Where T is in °K, the log Ap gives the decades of time that the 
relaxation curves should be shifted. The principle can be illustrated 
by referring to Figure A-1X where the test data corrected for density 
and temperature is shown on the left, and the log shift factor, log A^, 
is shown in the upper right. It can be seen from Figure A-ll that 
the log shift factor for -40° C is approximately +4, therefore the test 
data taken at -40°C should be shifted 4 decades to the right. The master 
curve in Figure A-ll does show that the test data taken at -40° C and 10° 
hours has been shifted to 10**^ hours. 


Another method (Reference 10) for determining the parameters for dynamic 
analysis is an emperical technique developed for determining G’ and G". 
The method involves vibrating a disk of the viscoelastic material and 
measuring its response. From this response, G 1 and G" can be mathe- 
matically determined. 

It was shown in Equations (43) and (44) that E 1 and E n are functions of 
frequency. Therefore, in a dynamic analysis the modulus used has to 
be based on the frequency of interest. In addition to frequency, the moduli 
also depend on strain, temperature, and pressure. The effects of each 
will now be discussed. 
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(Dynes/CM^) 







Effects of Fret 


The theoretical development of the dynamic moduli showed them to be 
functions of frequency. A master curve of dynamic shear moduli for a 
viscoelastic material is shown in Figure A~12. It is seen that the moduli 
increase with frequency at first and then decrease as the frequency is 
further increased. 

J} 

5.0| ’ 5. 


TOTE: Measurements made at 3 x 10**^ rms strain - -• 0 
levels. 

Log a™ = 0 at 25°C 


3.3.8 Effe 
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Figure A-13. Effects of Strain on Relaxation 
^Modulus 
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Figure A-14. Absolute Elastic Modulus vs. 

Amplitude of Sinusoidal Strain at 0. 5 cps 

3.3.9 Effects of Pressure 


An increase in pressure has also been shown to increase the dynamic 
moduli. Figure A-15 shows the increase in the relaxation modulus 
■with an increase in pressure from ambient to 1000 psi. The operational 
pressure of the SRB will be near 1000 psi. 
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Relaxation Modulus 
TP-H1011 Propellant 77°F 
-0. 06 + 0. 01 IN/IN 


Log t (MIN) 

Figure A- 15. Influence of Pressure on Relaxation Modulus 


3. 3. 10 Effects of Temperature 

A decrease in temperature causes an increase in the dynamic modulus. 
The increase in the relaxation modulus as temperature decreases is 
shown in Figure A-ll. Due to the rapid burning rate and the low thermal 
conductivity of the propellant, the temperature of the propellant will not 
increase significantly during its flight. Therefore, the temperature 
used for dynamic analysis should be the ambient temperature. 



4. 0 SUMMARY 


The physical properties of viscoelastic materials required for dynamic 
analysis has been defined and characterized. These properties are 
elastic modulus, shear modulus, Poisson’s ratio and the damping ratio. 

It was found that viscoelastic materials are incompressible and therefore 
Poisson’s ratio is 0. 5, It was also found that the dynamic elastic modulus 
and the dynamic shear modulus are complex and can be represented in the 
following form. 

E (-W) = E 1 (*>) + E" {“) ' (47) 

and G (-.«) = G 1 { B) + g" {« ) (48) 

It was shown that E l (w ) and E 11 (a) ) could be determined from the time 
response to a step strain, E r (t), or G'(u ) and G n (« ) could be determined 

by measuring the response of a vibrating disk. Dynamic shear modulus and 
dynamic elastic modulus are related by 


E(“) = 3G( U ) 

(49) 

The damping ratio, 5 , was found to be 

n 

6 = E 

(50) 

2E l 


The dependence of the storage moduli and the loss moduli on frequency, 
.. strain, .pressure, and temperature were discussed. __ All of these .effects 
must be considered during a dynamic analysis. 
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NASTRAN Computer Model 


The trapezoidal ring element is used to generate an axisymmetric model of a 
hollow cylinder to evaluate the thickness shear mode. The propellant is repre- 
sented by a 10 x 10 grid as shown in Figure B-l. The z dimensions of the ring 
elements are varied to obtain cylinders of different lengths. 

For data card preparation, only the cards which are different than those for 
normal mode analysis will be discussed. As mentioned previously, rigid format 
seven or ten can be used. In rigid format seven, the displacements of the grid points 
are used as the degrees of freedom, and in rigid format ten the normal vibration 
modes of the undamped system are used as the degrees of freedom. 


The METHOD card, which specifies an EIGR card in normal mode analysis, must 
be replaced by a CMETHOD card which specifies an EIGC card. The EIGC card 
selects the method for eigenvalue extraction and the search region. The specifi- 
cation of the search region was found to be the most critical area in complex 
eigenvalue analysis. If the search region is not specified with some degree of 
reasonableness, the computer time to extract the roots can become excessive or 
they cannot be found at all. The eigenvalue lies in the complex plane. Better 
results are obtained when a normal eigenvalue analysis is first made and an 
estimate is made of the decay rate based on the amount of damping. A typical 
search region is shown in Figure B-2. It is important to note that the search 
region for a system with less than critical damping lies in the second quadrant. 

The complex stiffness is calculated by NASTRAN by two different methods. All 
elements can be given equal damping by specifying G on a PARAM card or each 
element can be given a different value of damping by specifying GE on the element’s 
material property card. Viscous damping can be specified by using a scalar damper 
or a viscous rod element. 

One other property which must be input is Poisson’s ratio. As mentioned previously 
Poisson’s ratio for viscoelastic materials is 0. 5. NASTRAN cannot accept this 
value due to a 1-2 v term in the denominator of some equations. For this reason 
a value of 0, 495 is used. It is interesting to note that in the thickness shear mode 
motion only occurs in one direction, therefore, Poisson's ratio does not appear 
in the differential equations of motion. Computer runs were also made using a 
Poisson’s ratio of 0. 3. The eigenvalues did change by a small percentage probably 
due to radial motion still present in the computer model. 

A sample NASTRAN listing is given inAttachment. 
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FIGURE B-l. Computer Model 


ORIGINAL PAGE 1 8 
OP POOR QUALITY? 




T 


HIIMMMMH MMMMMMM mmm 

M'lHMMM'IMMMM'IMMUM/IMmiMMMMnMMimMMMMMrt 

MM M KM M M MMMMMMM M MM MMM MMMMMMMMMMMMMMMMMMMM / IBM 360-370 SERIES 

MMMMMMMKMMMMMMMMM m MMMMMMMMMMMMMMMMMMMMMMMMMMM // MODELS 91,95 

M M M M VI MM KMK"m M MMMM MMM MMMMM MMMM MMMM ‘■'MMMM MMMMMMMMMM /// 

MMM MMMMMM KMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM / /// 

"MMMMMMMMMMKMM ' M M M~M M MMMMM MM MM MM M M* MM MMM MM M MM M MMMM M M ///// 

MMMMMMMMMMMM MMHM MMMMMMMMMMMMMMMMMMHMMMMMMMMMMH ////// RIGID FORMAT ScRIES M 

MMMMKMMM M MMM“MMMMMM'MVMMMM MMMMMMMMMMMMMMM MMMMMfSMMMMMMMMM////// 

MMMMMM MMMMMM M M M MMM M M MM MMM MM MMMMMMMMMMMMMMMMMMMMMMMMMMMM////// 

mmmmmmmm'^mmmm 'wmmmmmmmmm'mmmmmmmm mmm mmmmmmmm m m*mm mm 7 MnM’ri M / / / / / / Lfc VLL 15. 1.0 

MMMMM MMMMMMM Mf- MMMMMM MMMMM MMMMMMM MM.MMMMKMMMM.MMMMMi’LMM//// / // / MM _____ 

v.MMMWmmMMMMMMM' MM'*MMM^MMMMMMM'fMMMMMMMM M ,M M M MMM MM K M MMM/// 7 ) '/ / / - M MM M 
MMMMM MMM MMMM V MMM MMMMMMM MMMM MM MMMMMM MMMMM MMMMMMMMM/ //////// --MMMMMMM 
' M.MMvmMMMMMMMMMMM MM MM MM MMMMM MMMMM MMM MMMM MM MMMM* ///7//7/// — MMMMMMMMM 
MM'-'MMMMMMMMM'-'MMM MMMMMMMM MMMHMMMMMM.M MMMM m ////// //mm- -mm.mmmm.mmmmm 

MMMMMMMM ' M v M7-‘M"m’ MMMMM M.MMMMMMMMM /MM — MM MM MM MM MMMMM MM MM MMMM 

__ MMM^MMMfM _MMMMMM MM MM MMM MMMMMM/// /// M MM — MMM MMM MMM M M MM MMMM MMMM _MM__ 

MMMMMVMMM 'MMMMM MMMMMM M 7/ / 7/7/7" MMMM MM — MMMMM MM'MMMM'm M MMM MM MMM MMMMM MM 

MMMMMMM U M M MMM MMMM// Ht /////MMM MMMM — MMMMMM MM MMM MM M _MMM M__ .MMM MM MMMM MM _ 

"mm "MM 777/ //7 MMMMMMMM ™'-m MMMMMMM 'MM MMM MM M MMM MM MMMM MM MMMM MM 

MM’-’M M' , 'M V M MMMM / /// ///MM MMMMMMMM M MMMMMM MMMMM MM MMMMMM M .MMMM MM MMMM MM 

•■MMMMMMM^" - ///// // M MMMMMM MMMM MMMMM MmTmMMMM M MMM MMMMMMMMMMM MM MMMMM 

MM/////// MMMMMM MMM MMMMMMM MMM M_MMMM MMMMMMM M MMM MM MM MM . MM MMMM 

* 7 / //MMMMMM M^MM MM MM MMMMM M " MM MMM MMMMM MM MM MMMM MMM'MMM MMMM MM 

., / / _ MM V MMMMMMMMMMM MM M MM MMMMMMM MMMMM M MMMMMMMHMMMMMMMM MHMMMMMMM MMMMMM MMM 

, 'mM w MM« v MM m MMMM'<MM MMMMMMM mmm MM-" MM MMMM MMMM MMMMMMMMM MMM'MMMMMMIIMMMMM 

^ MMMM.' M“MMMMMMmmM.‘1MMMMM mmmMH MMMMMMM.MMM MMMMMKMMMMM MM MMMMMM M M M MM MMMM 

“7mmv\ivxum\imm*Mm’m'mMMMM«MM MM MMMMMMMMMMM MMM MMMMMMMM M MM M MMMM M M MM M 

__MM MM MMMMMMMMMMM MMMMM KM — MMMMMM MMMMM MMMMMM MMMMMMMM. _MMMMM_ .MMMMMM 

MM MMMMMMM M V M MMMMMMM- ■MMMMMMM MMMMM M^MMMMMMMMMMMM" MMMMM MMMMMMM 

MMMMMMMMM MMMMMM — MMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM MMMMMMMMMMM 

K M M M MMMMMM MM u — -MMMMM y ' M M M M M MM m M M MM MMMMMMMMMMMMM~MM M K M KMM M MM 

_M y M MMMMMM MM^ MMMMMMM MMM M MMMMMM M M :1M MM M MMM MMM M MMMM MMMMMM SYSTEM GENERATION DATE - 5/13/72 

MMMMMMMM*-"- MmmMMMM"mmmm7mm7mmM~KMMMMMMm" ’ MMMM, MMMMMMMM 

MMM M M — MM MMMM ^ MM '■‘MMM, MMMMMMMM MMMMMMMMMMMMMKMMMKflM MMMMM 

" MM-- MM mmmMM MMM MMM" MMMMMMMMM M M KM MMMM KM MM M MMM KMK M M 

-MM MMMMMMM v MMMMMMM^MMMMMMMMMMMMMMKMMMMMMMMK.MMM 

MM7yKKK>.MMMMMHMKMMM 1 7KM,7M’KMKKKKMrHHHF7MK.MMMM 

- MMMMMMMMMMMMMM.MMMMMMMMMMMMMMMMMMMMM 

MMMMMMMM MMMMMMMMMMMMMMMMMMMMM 

MM MMMMM MMM MM MMM. MM 


££ 

S3 


APPENDIX -B ATTACHMENT' 


SEPTFMBCR 5> 1973 NASTRAN 5/13/72 


PAGt 


NASTRAN EXECUTIVE CONTROL DfcC'K ECHO 


ID SR8 LONGITUDINAL SHEAR ANALYSIS 

APP DISPLACEMENT 

T I «E 1 0 

SOL 7;<r 

CEND 








1 TITLE « NORMAL MODE ANALYSIS 

2 SPC#30 

3 CM ETHO O «L1 

-V ' 01 SP LACE ME NT a ALL 





NORMAL M'IDE ANALYSIS 


SCPJLMUfcR . 5, s L973^ „ NASTR AN \if 1J/ 72 


PAGt. 


3 


— 










CARD 



SORT 

f= 

6 B 

U L K 

DATA 

ECHO 


CfiUNT 

• 1 4 * 2 

« « 

3 .. 

4 

• • 

5 .. 

6 

7 •• 8 • i» 9 •• 10 • 



i- 

MRAPpi] ! 

1 

y 


, ii 

. 12 


..20 . 


2- 

CTRAPRG 2 

2 

3 


14 

U 


20 


3- 

CTUAPRO 3 

3 

4 


15 

14 


20 


4- 

CTPAPRG 4 

4 

5 


16 

15 


20 


5 — 

CTRAPRG 5 

5 

6 


17 

16 


20 


6- 

CTR4PRC. 6 

6 

7 


18 

17 


20 


7- 

CTRAPRG 7 

7 

a 


19 

18 


20 


8- 

CTRAPRG 8 

8 

9 


20 

19 


20 


Q — 

CTRAPRG 9 

9 

10 


21 

20 


20 


10- 

CTRAPRG 10 

10 

11 


22 

21 


20 


li- 

CTRAPRG 11 

12 



24 

23 


20 

§s 

12- 

CTPAPRG IZ 

13 

14 


25 

24 


20 


13- 

CTRAPRG 13 

14 

15 


26 

25 


20 


14- 

CTRAPRG 14 

15 

16 


27 

26 


20 

•Sj • 

- 

15- 

CTRAPRG 15 

16 

17 


28 

27 


20 

pig 

16- 

CTRAPRG 16 

17 

18 


29 

28 


20 

o ._ 

17- 

CTRAPRG 17 

18 

19 


30 

29 


20 


18- 

CTRAPRG 18 

19 

Ml 


31 

30 


20 

4g.,faf 

19- 

CTRAPRG 19 

20 

21 


32 

31 


20 


20- 

CTRAPRG 20 

21 

22 


33 

32 


20 


21- 

CTRAPRG 21 

23 

24 


35 

'34 


20 


2?- 

CTRAPPG 22 

24 

25 


36 

35 


20 

w 

23- 

CTPARRG 23 

25 

26 


37 

36 


20 

l 

24- 

CTRAPRG 24 

26 

27 


38 

37 


20 


2 5- 

CTRAPRG 25 

27 

28 


39 

38 


20 


26- 

CTPAPRG 26 

28 

29 


40 

39 


20 


27- 

CTRAPRG 27 

29 

30 


41 

40 


20 


28- 

CTRAPRG 28 

30 

31 


42 

41 


20 


29- 

CTRAPRG 29 

31 

32 


43 

42 


20 


30- 

CTRAPRG 30 

32 

33 


44 

43 


20 


31- 

CTRAPRG 31 

34 

35 


46 

• 45 


20 


32- 

CTRAPRG 32 

35 

36 


47 

46 


20 


33- 

CTPAPRG 33 

36 

37 


48 

47 


20 


34- 

CTPAPRG 34 

37 

38 


49 

40 


20 


35- 

CTRAPRG 35 

38 

39 


50 

49 


20 


36- 

CTRAPRG 36 

39 

40 


51 

50 


20 


3 7- 

CTPAPRG 37 

40 

41 


52 

51 


20 


38- 

CTRAPRG 38 

41 

42 


53 

52 


20 


39- 

CTRAPRG 39 

42 

43 


54 

53 


20 


40- 

CTRAPRG 40 

43 

44 


55 

54 


20 


41- 

CTRAPRG 41 

45 

46 


57 

56 


20 


42- 

CTPAPRG 42 

46 

47 


58 

57 


20 


43- 

CTPAPRG 43 

47 

4 8 


59 

58 


20 


44- 

CTRAPPG 44 

48 

49 


60 

59 


20 


45- 

CTRAPRG 45 

49 

50 


61 

60 


20 


46- 

CTRAPRG 46 

50 

51 


62 

61 


20 


47- 

CTRAPRG 47 

51 

52 


63 

62 


20 


48- 

CTPAPRG 48 

52 

53 


64 

63 


20 


49- 

CTRAPRG 49 

53 

54 


65 

64 


20 


50- 

CTPAPRG 50 

54 , 

55 


66 

65 


20 








SEPTE M BER 5, 1973 NASTRAN 5/13/72 


PAGt 


NORMAL MOPE AN ALYS IS 


4 





. S 

0 R T 

E 

D B U 

L K 

DATA 

ECHO 

CARD 









COUNT 

. 1 .. 2 

.. 3 

• * 

A 

.. 5 

• • 

6 •• 

7 .. 8 .. 9 .. 10 . 

51- 

CTRAPRG 51 

56 

57 


68 

67 


20 

52- 

CTRAPRG 52 

57 



69 

68 


20 

53- 

CTRAPRG 53 

58 



70 

69 


20 

5 A- 

CTPAPRG 5A 

59 

60 


71 

70 


20 

55- 

CTRAPRG 55 

60 

61 


72 

71 


20 

56- 

CTRAPRG 56 

61 

62 


73 

72 


20 

.O O 57- 

CTPAPRG 57 

62 

63 


7A 

73 


20 

53 - 

CTRAPRG 58 

63 

6A 


75 

7A 


20 

►0 S 59- 

CTPAPRG 59 

6 A 

65 


76 

75 


20 

O a 60- 

CTRAPRG 60 

65 

66 


77 

76 


20 

« 61- 

CTRAPRG 61 

67 

68 


79 

7 B 


20 

~TF“ ' 62- 

CTRAPRG 62 

68 

69 


eo 

79 


20 


CTRAPRG 63 

69 

70 


81 

80 


20 

6 A— 

CTPAPPG 6A 

70 

71 


82 

81 


20 

JHjQ 65- 

CTPAPRG 65 

71 

72 


83 

82 


20 

a fej 66- 

CTRAPRG 66 

72 

73 


8 A 

83 


20 

a *- 67- 

CTPAPRG 67 

73 

7A 


85 

8 A 


20 

42 «J 68- 

CT°APRG 68 

7A 

75 


86 

85 


20 

69- 

CTPAPRG 69 

75 

76 


87 

86 


20 

70- 

CTRAPRG 70 

76 



88 

87 


20 

71- 

CTRAPRG 71 

78 



90 

89 


20 

72- 

CTRAPRG 72 

79 

80 


91 

90 


20 

W 7 3- 

CTRAPRG 73 

80 

81 


92 

91 


20 

I 7A- 

CTPAPRG 7 A 

81 

82 


93 

92 


20 

03 75- 

CTRAPRG 75 

82 

83 


9A 

93 


20 

76- 

CTPAPRG 76 

83 

8A 


95 

9A 


20 

77- 

CTPAPRG 77 

8A 

85 


9 6 

95 


20 


CTRAPRG 78 

85 

86 


97 

96 


20 

79- 

CTRAPRG 79 

86 

87 


98 

' 97 


20 

80- 

CTRAPRG 80 

87 

88 


99 

. 98 


20 

8 1- 

CTRAPRG 81 

89 



101 

100 


20 

82- 

CTRAPRG 82 

90 

91 


102 

101 


20 

83- 

CTRAPRG 93 

91 

92 


103 

102 


20 

8A- 

CTPAPPG 8 A 

92 

93 


104 

103 


20 

85- 

CTRAPRG 85 

93 

9A 


105 

104 


20 

86- 

CTPAPRG 86 

9 A , 

95 


106 

105 


20 

87- 

CTRAPRG 87 

95 

96 


107 

106 


20 

83- 

CTRAPRG 98 

96 

97 


108 

107 


20 

89- 

CTRAPRG 89 

97 

98 


109 

108 


20 

90- 

CTRAPRG 90 

98 

99 


110 

109 


20 

91- 

CTRAPRG 91 

100 

101 


112 

1..U. 


20 .... 

92- 

CTRAPRG 92 

101 

102 


113 

112 


20 

93- 

CTRAPRG 93 


103 


1 1 A 

113 


20 

9A- 

CTRAPRG 9A 

103 

10A 


115 

i 1 A 


20 

95- 

CTPAPRG 95 

10 A 

105 


116 

115 


20 

96 - 

CTPAPPG 96 

105 

106 


117 

116 


20 

97- 

CTRAPRG 97 

106 

107 


1 IB 

117 


20 

93 - 

CTRAPRG 98 

107 

108 


119 

iia 


20 

99 - 

CTRAPRG 99 

106 

109 



120 

119 


20 


loo- ctraprg' ioo 109 no m "120 20 




NORMAL MODE ANALYSIS 


SEPTEMBER 


5, L973 NASTRAN 5/U/72 


PAOt 


SORTED BULK DATA ECHO 


CARD 


COUNT 

101- 

l 

nor. 

.. 2 . 

t 1 

• 3 . . 4 . . 5 . . 

INV MAX 

6 • * 

7 . . a 

• « 

9 .. 10 . 

ABC 1 


1 0 2— 

MIC 1 

-38.0 

lyo.o -38.0 24J.0 

20.0 

l 

1 

AUL2 


103- 

MIC 2 

-92.0 

500.0 -92.0 550.0 

20.0 

i 

l 

ABC3 


104- 

GBC3 

-104-0 

59070 -104.0 64670 

20.0 

l 

1 

ABC4 


10 5- 

SBC 4 

-338.0 

760.0 -138.0 810.0 

20.0 

i 

l 

ABC5 


106- 

EBC5 

-180.0 

995.0 -180.0 1045.0 

20.0 

l 

1 



107- 

GRDSET 




2456 




108- 

GRID 

l 

20. 0-0 0.0 






109- 

GRID 

2 

25. 0.0 0.0 






' L lO- 

GRID 

3 

30. 0.0 6.0 






ll 1- 

GRID 

4 

35. 0.0 0.0 






112- 

GRID 

5 

40. 0.0 0.0 




« 


113- 

GRID 

6 

45. 0.0 0.0 






lu- 

GP ID 

7 

50. 0.0 0.0 






ll 5- 

OR ID 

3 

55. 0.0 0.0 






lu- 

GRID 

9 

60. 0.0 0.0 






ll?- 

GRID 

10 

65. 0.0 0.0 






118- 

GRID 

11 

70. 0.0 0.0 

* 





119- 

GR ID 

12 

20. 0.0 30. 






120- 

GRID 

13 

25. 0.0 30. 






121- 

GR ID 

14 

30. 0.0 30. 






122- 

GRID 

15 

35. 0.0 3 0. 





w 

123- 

GRID 

16 

40. 0.0 30. 





1 

124- 

GR ID 

17 

45. 0-0 30. 




, 

vO 

125- 

GRID 

18 

50. 0.0 30. 




■ 



















PAGfc 


NOR HA L_HOD E ANALYSIS 


SEPTEMB ER 5,. 1973 N ASTRAN 5/13/7 2 



CARD 



S 

0 A T 

ED B 

ULK DATA ECHO 



COUNT 

. 1 

2 .. 

3 

« * 

4 .. 

5 .. 6 . . 7 .. 8 

.. 9 .. 10 . 


201- 

GRID 

94 


45. 

. 0.0 

240. 



202- 

GRID 

95 


50. 

0.0 

240. 



203- 

GR 10 

96 


55. 

0.0 

240. 



204- 

GRID 

97 


60. 

0.0 

240. 



205- 

GR ID 

°B 


65. 

0.0 

240. 



206- 

GRID 

99 


70. 

0.0 

240. 



207- 

GRID 

100 


20. 

0.0 

270. 


o o 

203- 

GR ID 

101 


25. 

o.o 

270. 



2 09- 

GRID 

102 


30. 

0.0 

270. 


"w'B” 

210- 

GRID 

103 


35. 

0-6 

270. 


© feJ. 

211- 

GRID 

104 


40. 

0.0 

270. 


2 S 

212- 

GRIO 

105 


45. 

0.0 

270. 

. 

W K 

213- 

GRID 

106 


50. 

0.0 

270. 

* 

£> tin 

214- 

grid 

107 


55. 

0.0 

270. 



215- 

GRIO 

108 


60. 

0.0 

270. 


El 

216- 

GR ID 

109 


65. 

0.0 

270. 



217- 

GRIO - 

110 


70. 

0.0 

270. 




... g,ff. 

210- 

GRID 

111 


20. 

o.o 

300. ‘ 


?|0r- 

r.» i o 

11? 



O.o 

100, 



"“220-“ 

'Grid 

113 


30. 

O.o’ 

300. 



221- 

GR ID 

114 


35. 

0.0 

300. 



222- 

GR ID 

115 


40. 

0.0 

300. 


w 

223- 

GRIO 

116 


45. 

0.0 

300. 


i 

224- 

GR ID 

117 


50. 

0.0 

300. 


H 

225- 

GRIO 

113 


55. 

0.0 

300. 



226- 

GR ID 

119 


60. 

o'.o 

300. 



227- 

GRID 

120 


65. 

0.0 

300. 



228- 

GRID 

121 


70. 

0.0 

300. 



229- 

MAT! 

20 


6244 

.0 0.495 1.658E-4 

0.364 


230- 

PARAM 

' COUPMASS1 







231- 

SPC1 

30 13 


11 

22 

33 44 55 

66 XYZ1 


232- 

GYZ1 

77 


88 

99 

110 121 



ENODATA 


♦ »NO ERRORS FOUND - EXECUTE NASTRAN PROGRAM** 








t'**USER 

INFORMATION 

MESSAGE 

302 8, 




38AR « 

23 C0AR » 

OR# 

46 








PAGE 


NORMAL MOD E_ A N A L YS_I_S 


SEPTEMBER 5, 1973 NAS TRA N 5/13/72 


a 


COMPLEX EIGENVALUE ANALYSIS SUMMARY ^INVERSE POWER WITH SH1FTS< 




NUMBER OF EIGENVALUES EXTRACTED 

1 





NUMBER OF STARTING POINTS USED 

1 




o o 

NUMRFfc OF STARTING POINT OR SHIFT POINT MOVES „ . 

0 




- Si 

TOTAL NUMBER OF TRIANGULAR DECOMPOSITIONS .... 

1 





TOTAL NUMBER OF VECTOR ITERATIONS 

10 






ROOT 

NO. 

EXTRACTION 

ORDER 

EIGENVALUE 

*REAl< ? I MAG< 


FREQUENCY 

%CYCLES< 

damping 

COEFFICIENT 

1 

1 

-4-2993656 01 

2.4381Q4E 

02 

3.8603626 01 

3.52b81uL-0*. 

2 

2 

-8.9612006 01 

5.0Q1758E 

02 

8 . 087868 E 01 

3. 52ballfc-01 

3 

3 

- 1.04612 IE 02 

5.932393b 

02 

9.4416966 01 

3.5268106-01 

4 

4 

-1.4506096 02 

6. 226182E 

02 

1.3092386 02 

3- 526olOt— 01 

5 

5 

-1.800453E 02 

1.021009E 

03 

1. 624987b 02 

3. 526 81Ui.— 01 





NORMAL MODE ANALYSIS 


sfptemrep 


5, 1973 NASTRAN S/li/72 


fMiu 


11 ) 


COMPLEX EIGENVALUE 

« -4.299365E 01, 

2.438104E 

02 





* 


C 

0 M P L 6 X 

EIGENVECTOR NO. 

*REAL/IMAGlNAfiY< 


1 


POINT 10. 

TYPE 

T l 

T2 

T3 

R1 


R2 

R3 

L 

G 

-1.218973T-01 

0.0 

9.805287E-01 

0.0 

0.0 


0.0 



6. 4696? IE-07 

0.0 

-3.207562E-07 

0.0 

0.0 


0.0 

2 

G 

-1.008059E-01 
4.763 742E-07 

0.0 

0.0 

9.631643E-01 

-5.12304BE-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

3 

r. 

-8. 5558536-02 

0.0 

9.144636E-01 

0.0 

0.0 


0. 0 



3 , 46002 26-0 7 


— 5. 6827 92 E -07 

0.0 

0.0 


0.0 

4 

G 

-7.362592E-02 
2. 475405C-07 

0.0 

0.0 

8.434793E-01 
— 5.456513E— 07 

0.0 

0.0 

0.0 

0.0 

< 

0.0 

0.0 

5 

G 

-6.309706E-02 

0.0 

7. 545497C-01 

0.0 

0.0 


0.0 



1. 704041E-07 

0.0 

-4.6965356-07 

0.0 

0.0 


0.0 

6 

G 

-5.2829 1 7E-02 
1. 107600E-07 

0.0 

0.0 

6.511883E-01 

-3.607460E-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

7 

G 

-4. 2179583-02 

0.0 

5.365742E-0 1 

0.0 

0.0 


0.0 

W 


6. 585 1 1 7E-08 

0.0 

-2.386519E-07 

0.0 

0.0 


0.0 

J-4 ^ 

i_rt 

G 

-3.08975 8 F-02 
3.38624SE-08 

0.0 

0.0 

, 4.13535Bfc-01 
-1.229655E-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

9 

G 

-1.915 107E-02 

0.0 

2. 841575E-01 

0.0 

0.0 


0.0 



1 . 3342 67E-03 

0.0 

-3.272350E-08 

0.0 

0.0 


0.0 

10 

G 

-7. 7358 LtC-03 
2 . 77241 9E-09 

0.0 

0.0 

1.488124E-0I 

1.406336E-08 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

11 

G 

0.0 

0.0 

o.-o 

0.0 

0.0 


0.0 



0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

i"2 

G 

-2. 562813E— 02 - 
8. 32646 5F- 07 

0.0 

0.0 

l.OOOOOOE 00 
0.0 

0.0 

0.0 

0.0 

0.0 


0.0 
0.0 . 

13 

G 

-2.1303 7AF-02 

0.0 

9. 709264E-0 1 

0.0 

0.0 


0.0 



6. 00 7 166b - 07 

0.0 

— l .073025L-07 

0.0 

0.0 


0.0 

""14 

G 

-1.81551 5E-02 
4. 3110676-07 

0.0 

0.0 

9. 174005E-01 
-1.630101E-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

15 

G 

-1.55tl2B=-02 

0.0 

8.405590L-01 

0.0 

0.0 


0.0 



3.0170416-07 

0.0 

—1.6627536—07 

0.0 

0-0 


0.0 

16 

G 

-1.293326E-02 
2. 0186676-07 

0.0 

0.0 

7.444257E-01 

-1.3203316-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 









N09M4L MODE ANALYSIS 


SEPTEMBER 5, 1973 NA STRAN 5/13/72 PAGE 1 1 


COMPLEX eigenvalue 

ft -4.29 9365 E 01, 

2.438104E 

02 





* 



C 

0 M P L E X 

eigenvector no. 

*REAL/1MAGINARY< 


1 


POINT 

ID. 

TYPE 

T 1 

T2 

T3 

R1 


R2 

R3 


17 

G 

-1.03386IE-02 

0.0 

6.330248E-01 

0.0 

0.0 


0.0 • 




1. 2 577 8 7C-07 

0.0 

-7.5778306-08 

0-0 

0.0 


0.0 


18 

G 

— 7 . 50 1 74 7E-03 
fc. 99232 6C-08 

0.0 

0.0 

5. 10418 7C-01 
-1.3515B3E-08 

0.0 

0.0 

IIK&til 


0.0 

0.0 

O O 
*9 # 

19 

G 

-4.5606206-03 

0.0 

3.807130E-01 

0.0 

0.0 


0.0 

% b 
© e 



3.177 63 IE-08 

0.0 

3.886636E-08 

0.0 

0.0 


0.0 


20 

G 

-1.8 1206 3E-03 
9. 097 075C-09 

0.0 
0 . 0 

2.48 1694E-0 1 
6.649805E-00 

0.0 

0.0 

0.0 

0.0 


. 0.0 
0.0 

*TJ 

5? hi 

21 

G 

7.491 151F-05 

0.0 

1. 174533E-01 

0.0 

0.0 


0.0 




-5. 1524056-10 

0.0 • 

5.659170E-08 

0.0 

0.0 


0.0 

Sgg 

?2 

0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o-o 

0.0 

0.0 

0.0 


0.0 

0.0 


23 

G 

-3. 1451 44E-03 

0.0 

9. 9232 75E-01 

0.0 

0.0 


0.0 

w 



1.268 146E-06 

0.0 

6.432G53E-07 

0.0 

0.0 


0.0 

~r — 

o 

24 

G 

- 1. 049 323L-03 
9.071 106E-07 

0.0 

0.0 

9.6874386-01 

5.256587C-07 

0-0 

0.0 

0.0 

0.0 


0.0 

0.0 


25 

G 

-3. 896000E-05 

0.0 

9. 149157E-01 

0.0 

0.0 


0.0 




6. 4914716-07 

0.0 

4. 5291556-07 

0.0 

0.0 


0.0 


26 

G 

5 .926276 r -Q4 
4.5433436-07 

0.0 

0.0 

6.372628E-01 

4.047499C-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


27 

G 

1.0 19964C-03 

0.0 

7 . 399346E-01 

0.0 

0.0 


0.0 




3.04695 IE-07 ' 

0.0 

3.733484E-07 

,0.0 

0.0 


0.0 


28 

G 

1.2918006-03 
1 .908 6306-07 

0.0 

0.0 

6.273241E-01 
3. 489028E-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


29 

G 

L.398P94E-03 

0.0 

5.039493E-01 

0.0 

o.o 


0.0 




1.07221 96-07 

0.0 

3.21 13 70E-07 

0.0 

0.0 


0.0 


30 

G 

1.289513E-03 
4. 976 563E-08 

0.0 

0.0 

3. 7453286-01 
2. 8002 3 3 E— 07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


3 l. 

G 

8 . 899 50 2E-04 

0.0 

2.4405906-01 

0.0 

0.0 


0.0 




1. 504150E-03 

0.0 

2.167030E-07 

0-0 

0.0 


0.0 


32 

G 

2.1 98 S74E- 04 
-3.7033916-10 

0 . 0 
0.0 

1. 164155E-01 
1.2448416-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 
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COMPLEX EIGENVALU E tf -4.2993656 01 , 2 .438 1 04E 02 

COMPLEX EIGENVECTOR NO. 1 

?REAL /IMAGINARY < 


POINT ID. 

TYPE 

T 1 

T 2 

T3 

Rl 

R2 

R3 


33 

G 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 




0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


3A 

G 

I . 1 49066E-03 
I. 682043, = -06 

0.0 

0.0 

9.794863E-01 

1.554316E-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


1 

1 

j 

! 

w 

1 

G 

2. 395362E-03 

0.0 

9. 5979 86E-01 

0.0 

0.0 

0.0 


• 


1.2057035-06 

0.0 

1.436731E-06 

0.0 

0.0 

0.0 


36 

G 

3. ll45S°F-03 
8.630S89E-07 

0.0 

0.0 

9. 090770E-01 
1.324316E-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


_37 

G 

3. 3858605-03 

0.0 

8.333428E-01 

0.0 

0.0 

0.0 




6.0 46 26 25 -07 

0.0 

1.209190E-06 

0-0 

0.0 

0.0 


38 

G 

3.330474E-03 

4.0625835-07 

0.0 

0.0 

7 .377 120E-01 
1.08 43 50 E-06 


0.0 

0.0 

0.0 

0.0 


39 


2.999298E-03 

0.0 

6.267717E-01 

0.0 

o’.o 

0.0 


W 


2.5 5?0 756-07 

0.0 

9.471724E-07 

0.0 

0.0 

0.0 


^ AO 

G 

2. 4207525-03 
1.4417555-07 

0.0 

0.0 

5.05I230E-01 

7.947182E-07 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


41 

G 

1. 675395C-03 

0.0 

3.7740776-01 

0.0 

0.0 

0.0 




6. 749 5 7 65- 00 

0.0 

6.245563C-07 

0.0 

0.0 

0.0 



G 

8.6156355-04 

2.0765255-08 

0.0 

0.0 

2.482598E-01 
4.3 52262E-07 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


_ A3 

G 

2.6021555-04 

0.0 

1.2138 72E-01 

0.0 

0.0 

0.0 




-3.0233065-10 

0.0 

2.266545E-07 

0.0 

0.0 

0.0 


44 

G 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


45 

G 

8. 181 62 OE-04 

0.0 

9.708875E-0I 

0.0 

0.0 

0.0 




1^9635255-06 

0.0 

2. 6547356-06 

0.6 

0.0 

0.0 


46 

G 

1.4491635-03 

1.4094115-06 

0.0 

0.0 

9.525841E-01 

2.556211E-06 

0.0 

0.0 

0.0 

0.0 

0.0“ 

0.0 


47 

G 

1.8582145-03 

0.0 

9.039427E-0 1 

0.0 

0.0 

0.0 




1.0096935-06 

0.0 

2^4045146-06 

0.0 

0.0 

0.0 


48 

G 

2.02622 5E-03 
7. 073825E— 07 

0.0 

0.0 

8. 3038 79E-01 
2.205611E-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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COMPLEX EIGENVALUE 

H -4.299365C 01, 

2.438104E 

02 








C 

0 M P L E X 

EIGENVECTOR NO. 

*REAL/IMAGINARY< 


1 



POINT ID. 

TYPE 

T 1 

T2 

T3 

R 1 


R2 

R3 


49 

G 

1.962661S-03 

0.0 

7 .367498 E-01 

0.0 

0.0 


0.0 




4. 7526766-07 

0.0 

1.965017E-06 

0.0 

0.0 


0.0 


50 

,0-0 — - 

G 

1. 706339E-03 
?.9fu,‘mr~07 

0.0 
0. 0 

6.2763296-01 

1.6niV,ll?C”06 

0.0 

0.0 

0.0 

0.0 


o.o 

o.() 


*3 W 

J*-g3 - 51 — 

G 

1 .3060096—03 

0.0 

5. 07 49 8 8 E-01 

0.0 

0.0 


0.0 


§|_ 


1. 6374055-07 

0.0 

1.3S2300E-06 

0.0 

0.0 


0.0 


52 

G 

8.366492S-04 
7. 906692E-08 

0.0 

0.0 

3. 807095E— 01 
1.0531 10E-06 

0.0 

0.0 



0.0 

0.0 


.l£__ 53 

G 

4. 0082°4S-04 

0.0 

2. 513349E-01 

0.0 

0.0 


0.0 




2 .4382566-08 

0.0 

7 . 0 7946CE-07 

0.0 

0.0 


0.0 


|^*573ESoS2 

G 

1.0759366-04 

-2.843490E-10 

0.0 

0.0 

1.232286C-01 
3. 542937E-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


55 

G 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

tri 



0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

-OX 

56 

G 

-3.939645E-06 
2 . 06044 56-06 

0.0 

0.0 

9. 6801406-01 
3 . 864448E-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


57 

G 

-2.825185C-06 

0.0 

9. 5003 89E-01 

0.0 

0.0 


0.0 




1.4775835-06 

0.0 

3.792674E-06 

0.0 

0.0 


0.0 


58 

G 

-2.020747C-06 

.1.0568636-06 

0.0 

0.0 

9.0204685-01 
3.60 10 53E-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


59 

G 

_=i.A124A5 6-06 0.0 

8. 29 30 8 3 E-01 

0.0 

0..0 


0.0 




7 . 38741 6F- 07 

0.0 

3.310635E-06 

0.0 

0.0 


0.0 


60 

G 

-9.4616295-07 

4.948533S-07 

0.0 

0.0 

7.364896E-01 
2 .9400 58 £-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


61 

G 

-5.925094C-07 

0.0 

6. 2808366-01 

0.0 

0,i) 


0.0 




3.0989066-07 

0.0 

2. 507265E-06 

0.0 

0.0 


0.0 


62 

G 

-3. 3329336-07 
1.7431845-07 

0.0 

0.0 

5 . 08 4775E-0 1 
2. 029777E-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 


63 

G 

-1. 552484E-07 

0.0 

3.8 18755E-01 

0.0 

o.o 


0.0 




8. 1 198725-08 

0.0 

1. 524380E-O6 

0.0 

0.0 


0.0 


64 

G 

-4. 7282106-08 
2.4730506-08 

0.0 

0.0 

2. 523409E-01 
1.007290E-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 
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COMPLEX FI G5NVALUE 

« -4.2993655 01, 

2.4381045 

02 







C 

0 M PL E X 

EIGENVECTOR NO. 

2!RE AL/ I MAG I NARY < 


1 


I'OINI It), 

type 

Tl 

n 

f 3 

in 


n 

R3 

65 

G 

l. 0390095-09 

0.0 

1.237865E-01 

0.0 

0.0 


0.0 



-S.427849E-10 

0.0 

4. 941268E-07 

0.0 

0.0 


0.0 

66 

g" 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

67 

G 

-8.2565625-04 

0.0 

9. 7080285-01 

0.0 

0.0 


0.0 



1.956982c-06 

0.0 

5.097098 1-06 

0.0 

0.0 


0.0 

60 

g 

-1^4545305-03 
1. 3978325-06 

0.0 

0.0 

9.525793E— 01 
5.O49449L-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

69 

G 

-1.8620475-03 

0.0 

9. 039381 E-Ol 

0.0 

0.0 


0.0 * 



9.948490E-07 

0.0 

4.8 127 19E-06 

0.0 

0.0 


0.0 

70 

G 

-2. 028900E-03 
6.9119935-07 

0.0 

0.0 

8.303837E-01 

4.424269E-06 

0.0' 
0 . 0 

0.0 

0.0 


0.0 

0.0 

71 

G 

—1.9644495-03 

0.0 

7.367461E-01 

0.0 

0.0 


0.0 

W 


4. 595932E-07 

0.0 

3.917166E-06 

0.0 

0.0 


0.0 

^ 72 

o 

G 

-1.7074655-03 

2.8506675-07 

0.0 

0.0 

6.276297E-01 

3.3224415-06 

0.0 

0.0 

0.0 

0.0 


0.0 
0 . 0 

73 

G 

-1.3067145-03 

0.0 

5.0749645-01 

0.0 

0.0 


0.0 

• 


1 . 533 1 02E-07 

0.0 

2.669432E-06 

0.0 

0.0 


0.0 

74 

G 

-3. 3693S9L-04 
7.2385385-03 

0.0 

0.0 

3.807077E-01 

1.986332E-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

75 

G 

-4.0891635-04 

0.0 

2 . 5 1333 SE-0 1 

0.0 

0-0 


0.0 



2.11 1750c- 03 

0.0 

1. 298597L-06 

0.0 

0.0 


0.0 

76 

G 

-1.0759035-04 
-1 . 143668F-09 

0.0 

0.0 

1.232281E-01 
6 . 29 5029C-07 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

77 

G 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 



0.0 

0.0 

0.0 

0.0 

0.0 

• 

0.0 

78 

G 

-1. 155431E-03 
1 .6729055-06 

0.0 

0.0 

9.7947735-01 

6.2661795-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

79 

G 

-2.3999365-03 

0.0 

9.597895E-01 

0.0 

0.0 


0.0 



1. 1Q6606S-06 

0.0 

6.226535E-06 

0.0 

0.0 


0.0 

80 

G 

-3 . 117 8 4 IE-03 
8. 3324355-07 

0.0 

0.0 

9 . 09068 2E-0 1 
5.9334135-06 

r 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 
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COMPLEX . EI GEN VALUE » -4.2993656 01, 2.438104E 02 


COMPLEX EIGENVECTOR NO* 1 

?R E AL/IMAG I NARY < 



POINT 10. 

TYPE 

T 1 

T2 

T3 

PI 

R2 

R3 


81 

G 

-3. 39812 IE-03 

0.0 

8 . 33 3347E-0 1 

0.0 

0.0 

0.0 




5. 7761043-07 

0.0 

5.4442863-06 

0.0 

0.0 

0.0 


82 

G 

-3.33I977E-03 
3. 796906F-07 

0.0 

0.0 

7. 377049E-01 
4.805522E-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

O < 

9 ' 88 

G 

-3.000228E-03 

0.0 

6.2676586-01 

0.0 

0.0 

0.0 

~ r ^J ! 

Sj 

g 

£J5 


2. 3 13700E-07 

0.0 

4.0568813-06 

0.0 

0.0 

0.0 


G 

-2.4292663-03 
1 . 2478983-07 

0.0 

0.0 

5. 05 118 3 £-01 
3. 23Q052E-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

JSJ 

a 

85 

G 

-1. 675627E-03 

0.0 

3.7740436-0 1 

0.0 

0.0 

0.0* 

L 


5.4122103-08 

0.0 

2.3885293-06 

0.0 

0.0 

0.0 

a *w 86 

G 

-8 . 6 16298E-04 
l. 3887433-09 

0.0 

0.0 

2.482576E-01 
l .546768 E-06 

0.0 ’ 
0.0 

0.0 

0.0 

0.0 

0.0 


87 

. 0 


0.0 . 

_ ).?n06300l 

0.0 

0.0 

0.0 

w 


-2. 3 19 7 893-09 

0.0 

7.4243973-07 

0.0 

0.0 

0.0 

ro 88 

o 

G 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


a9 

G 

3. 140246E-03 

0.0 

_._ 9. 9231483-01 

0.0 

0.0 

0.0 




1.293364E-06 

0.0 

7.279742E-06 

0.0 

0.0 

o.-o 


90 

G 

1. 045843E-03 
9.1 55754F-07 

0.0 

0.0 

9.687310E-01 

7.2090346-06 

0.0 

0.0 

, 0.0 
0.0 

0.0 

0.0 


91 

G 

• 3.6476791-05 

0.0 

9.1490353-01 0.0 

0.0 

0.0 




6. 495297E-07 

0.0 

6.851943E-06 

0.0 

0.0 

0.0 


92 

G 

-5.9435613-04 

4.496609E-07 

0.0 

0.0 

8.372515E-01 

6.2800453-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


93 

G 

-1. 021 1133-03 

0.0 

7.3992483-01 

0.0 

0.0 

0.0 




2.9659746-07 

0.0 

5.5342963-06 

0.0 

0.0 

0.0 


96 

G 

-1.292510E-03 
1 . 805 830F-07 

0.0 

0.0 

6.273159E-01 

4.6595943-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


95 

G 

— 1 . 3992 83 c -03 

0.0 

5. 03942 8E-01 

.0.0 

0.0 

0.0 




9.6075726-09 

0.0 

3.702295E-06 

0.0 

0.0 

0.0 


96 

G 

-1.2896833-03 

3.9483286-08 

0.0 

0.0 

3. 745281E-01 
2.7101363-06 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 
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COMPLEX EIGENVALUE 

« -4.299365c 01, 

2.438104E 02 










C 

0 M P L E X E 

IGENVEOTOR NO. 

3REAL/ 1 MAGINARY< 


1 




POINT 10. 

TYPE 

T 1 

T2 

T3 

R1 


R2 

R3 



97 

G 

-8.899940E-04 

0.0 

2 .44056 IE-0 l 

0.0 

0.0 


0.0 





7.942024E-09 

o.o 

1.731771E-06 

o.o 

0.0 


0.0 



98 

G 

-2. 198 826E-04 
-2. 12442 IE-09 

0.0 

0.0 

1.184142E-01 

8.206930E-07 

0.0 

o.o 

0.0 

0.0 


0.0 

0.0 



99 

G 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 





0.0 

0.0 ' 

0.0 

0.0 

0.0 


0.0 



100 

G 

2.562455E-02 
1 .037 440E-06 

0.0 

0.0 

9.999047E-01 

7.984262E-06 

0-0 

0.0 

0.0 

0.0 


o.o 

0.0 



101 

r. 

non it-o? 

0.0 

9.7091 12C-01 

o.o 

0.0 


0.0 





7.7096 U1E-07 

0. 0 

7 . 8594 L 511-06 

0.0 

0.0 


0.0 



102 

G 

L.815323E-02 
5.7618 l 5E-07 

. 0.0 
0.0 

9. L73859E-01 
7.487722t-06 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 



103 

G 

1. 550989E-02 

0.0 

8.405455l-01 

0.0 

0.0 


0.0 


w 



4. 2564865-07 

0.0 

6.877421E-06 

0.0 

0.0 


0.0 


r 

to 

v— ■ 

104 

G 

1.298229E-02 

3.056070E-07 

0.0 

0.0 

7.444139E-01 

6.075591E-06 

0.0 

0.0 

0.0 

0.0 


o.o 

0.0 



105 

G 

1.033797E-02 

0.0 

6. 330148E-01 

0.0 

0.0 


0.0 





2.083842 c -07 

0.0 

5.129852E-06 

0.0 

0.0 


0.0 

■ 


106 

r, 

7.501367E-03 
1.298604E— 07 

0.0 

0.0 

5. 104108t“01 
4.088665E-06 

0.0 

0.0 

0.0 

0.0 


o.o 

0.0 



107 

G 

4. 56042 6~ -03 

0.0 

3-807074E-01 

0.0 

0.0 


0.0 





6.821 460E- 08 

0.0 

3.000698C-06 

0.0 

o.o 


0.0 



108 

G 

1. 3 12001E-03 
2.3577285-03 

o.o 

0.0 
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APPENDIX C 

SRB SEGMENT TEST REQUIREMENTS 



APPENDIX C 


TITLE: SRB Segment Test Requirements 
1. 0 PURPOSE 


The purpose of conducting tests on segments of the SRB is to evaluate 
the analysis of prcpellant modes. . 

2.0 OBJECTIVE 


The objective of the segment tests is to determine the resonant frequencies, 
deflection responses, mode shapes and stress distribution in the first two 
longitudinal axisyrnrnetric thickness shear modes at several burn times, 

3.0 TEST SPECIMEN REQUIREMENTS 


The test segment shall be fabricated from a typical motor segment case. 

The dynamic properties ’of the propellant shall be identical to the SRB 
operational propellant. Grain design including liner and flap design shall 
be representative of actual SRB flight hardware at a forward segment 
representative of the smallest bore diameter. 

The length of the segment shall be approximately 300 inch. If SRB 
operational segment length is less than 300 inch, it shall determine the 
length of the test segment. 

The propellant shall be configured to several different bore diameters repre 
sentative of motor segment conditions or burn times as follows: 

1) Preignition 

2) Fundamental acoustic combustion resonance equal to fundamental 
longitudinal shear 

3) A bore diameter which yields the highest strains (See Section 4. 4. 3) 

4. 0 TEST ENVIRONMENT 


The vibratory excitation shall be sinusoidal sweep. The test frequency 
range shall be 10 to 50 Hertz. The sweep rate shall be one octave per 
minute. 


C-l 



Vibratory excitation forces shall be applied longitudinally to the motor 
segment case at the fore and aft interface. The forces shall be uni- 
formity distributed at the case/exciter interface. The net force at each 
interface shall be equal in magnitude, in phase and sufficient to cause 
(1) a Ig peak sinusoidal acceleration and (2) an excitation equal to the 
maximum expected levels at the interfaces for each burn time. 

The test segment shall be maintained at a uniform ambient temperature 
for a sufficient time prior to starting test so that a uniform temperature 
through-out the segment is obtained. This same ambient temperature 
shall be maintained during the test. 

5.0 INSTRUMENTATION 


Test input and response vibration levels shall be measured with accelero- 
meters. A sufficient number of accelerometers shall be placed on the 
fore and aft motor case interfaces to monitor the motion. Accelerometers 
shall also be attached to the surface and embedded within the propellant 
to evaluate the fundamental mode shape through-out a cross-section formed 
by the longitudinal axis and radius vector. Measurements shall be made at 
grid points formed by radii and longitudinal lines within the plane. . 

Strain gage measurements shall be made at the same grid point measure- 
ments defined for accelerometers so that the strain distribution can be 
evaluated. 

6. 0 DATE - REQUIREMENTS 

Data shall be taken during each test to insure that frequency response, 
mode shapes, acceleration and strain levels can be evaluated. 
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References and Bibliography 


The bibliography consists of two sections. Section 1 contains those 
references which were found to be very useful. These references are 
abstracted for future reference. Section 2 contains those references which 
were used, but not considered as applicable as the references in Section 1. 

Section 1. 

Achenbach, J. D. : "Axial Shear Vibrations of a Cylinder of 
Decreasing Thickness, 11 AIAA Journal, Vol. 4, No. 7, pp. 1233-1240, 

July 1966. 

This paper presents an analytical study of the axial shear vibrations 
of a long hollow cylinder that is subjected to time-dependent body 
forces in the axial direction. The outer surface of the cylinder is 
bonded to a rigid case, and the inner radius increases monotonically 
with time. An expression is determined for the shear stress at the 
bond interface. It is shown that the instantaneous frequency of the 
shear bond stress increases, and that its amplitude decreases toward 
burnout time. The axial shear vibrations of an ablating viscoelastic 
cylinder are discussed briefly . 

Anon. : "Final Report Frequency Investigation of the TU-122-1511. 11 
Minuteman Motor, " TR Thiokol Chemical Corporation. 

This report investigates the input to the first stage Minuteman motor 
during transportation. The critical speed was determined for traversing 
road course with different height boards. 

Anon. : "Final Report Resonance Search and Response Test Stage 1 
Minuteman Motor with Embedded Instrumentation (TU-122-1834. 1062, TTM 009), 
TR Thiokol Chemical Corporation. 

^•Anon: "Final Report Resonance Search % Stage 1 Minuteman Motor 
(TU-122-1834. 633; TTM-010)" TR Thiokol Chemical Corporation. 

These test reports identified the resonant frequencies of the first stage 
Minuteman motor which could be excited during transportation and 
handling. The frequencies identified were the coupled thickness shear 
mode, the lobar mode, and lateral bending modes. Also identified 
were the resonant frequencies cf the dome and nozzle. Responses were 
measured at the various resonances with embedded instrumentation. 
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^‘Anon. : "Ground Dynamic Tests (WS-133A)", Test No. T2-1704, 
Boeing Airplane Company, Seattle, Washington, January 1961. 

Tests were conducted on a full-scale Minuteman missile, containing 
inert propellant, to determine structural dynamic characteristics 
of the airborne vehicle. 

Objectives of the Test were* 

1. To define modal characteristics (frequencies, mode shapes 
and damping); 

i 

2. To obtain slope measurements at anticipated flight control 
sensor locations; 

3. To monitor auxiliary equipment resonances; 

4. To measure the missile’s response to a lateral thrust component, 
simulated by a sinusoidal force applied at the nozzle hinge. 

Baltrukonis, J. H. and Gottenberg, W. G. ; "Thickness -Shear 
Vibrations of Circular Bars, " Space Technology Laboratories, Inc. Report 
No. GM-TR-0165-00518, Nov. 1958. 

Exact general solutions of the three-dimensional elasticity 
equations of motion in polar cylindrical coordinates are obtained 
for axisymmetric and antisymmetric thickness-shear vibrations. 
These solutions are applied in solving five solid and hollow circular 
bar problems. Natural frequencies are tabulated and mode shapes 
are plotted. 

^ ' Baltrukonis, J. H. : "A Survey of Structural Dynamics of Solid 
Propellant Rocket Motors, " NASA CR-658, Dec. 1966. 

In the paper a survey of the dynamic problems of solid propellant 
rocket motors is presented starting from the simplest model thereof 
and proceeding, step-by-step, to the consideration of more sophisti- 
cated and realistic models. The consideration is restricted to 
infinitesimal deformations of propellant grains with linear mechanical 
properties. Substantial progress has been achieved towards the 
solution of many important dynamical problems. Pertinent develop- 
ments are summarized and indications are made along which lines 
future study shouM proceed. A few illustrative solutions are 
included in areas wherein progress has been substantial. 
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8* Baltrukonis, J. H. j Gottenberg, ~W. G. ; and Schreiner, R. N. : 

"Axial-Shear Vibrations of an Infinitely Long Composite Circular Cylinder," 

J, Acoust. Soc. AM. , Vol. 33, No. 11, Nov. 1961, pp. 1447-1457. 

Exact general solutions of the three-dimensional elasticity equations 
of motion in polar cylindrical coordinates are written for axisymmetric 
axial-shear vibrations. The frequency equation follows immediately 
from the boundary conditions for the problem of the infinitely long, 
composite cylinder with two concentric circular -cylindrical layers 
which are perfectly bonded at their interface. The branches of the 
frequency equation are plotted and analyzed. The conditions are 
pointed out under which it is possible to obtain reasonably accurate 
estimates of the natural frequencies by assuming that the motions 
of the casing and core are uncoupled. 

Beyer, E. B. : "Nonlinear Mechanical Behavior of Solid Propellants, " 

AIAA Paper No. 65-159. 

The mechanical response of standard solid propellants was studied 
by extensive analysis of data measured under conditions of constant 
strain rate, constant strain, and dynamic shear strain. Nonlinear 
viscoelasticity has been found to occur when propellant samples are 
strained beyond a few tenths of one percent by tensile test methods 
currently used by most investigators. Studies conducted over a wide 
range of strain rates {10 ^ to 10 ra in “ ~ ) indicate that nonlinearity can 
occur (1) by loss of reinforcement cue to dewetting and (2) by the 
"Mullins Effect" in a matrix with chemical adhesive bonding between 
binder and filler. In case (1) dewetring was observed to depend only 
on the applied stress and the temperature. The linear viscoelastic 
response obtained from small consrant strain rate and dynamic data 
differed from the constant strain (2 ro 5 percent) stress-relaxation 
modulus by as much as an order of magnitude. At very low strain rates 
equilibrium behavior was obtained cp to strains as high as 7 percent. 
The time and temperature dependence of both the reinforced and non- 
reinforced modulus is discussed and related to long term ambient tests 
and actual motor behavior. 

10 

* CPIA Publication No. 21: "Solid. Propellant Mechanical Behavior 

Manual/* 

This document is a general information reference on mechanical 
properties of solid propellants. It covers mechanics of viscoelasticity, 
failure criteria, procedures for physical characterization, strain 
measurement techniques nomenclature, symbology, surveillance for 
reliability and aging, crosslink density and equilibrium stress-strain 
measurements. 
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' CpiA Publication No. 214: "Handbook for the Engineering 
Structural Analysis of Solid Propellants, " Prepared by Fitzgerald, J. E. 
andHufferd, W. L., May 1971. 

This Handbook represents an attempt to present an accurate report 
and evaluation of the current state-of-the-art of solid propellant 
grain structural integrity. Major emphasis is given to the require- 
ments for meaningful material characterization, structural analysis 
and failure analysis of solid propellants. The Handbook has been 
written primarily for the designer/analyst. 

T ? 

• Davey, A. B. and Payne, A. R. : Rubber in Engineering Practice , 
PaLmerton Publishing Co. , Inc. , New York, 1964. 

This text describes the behavior of rubber under rapidly changing 
stresses such as vibrations and shocks, and of the dependence of 
this behavior on such factors as temperature, frequency and rate 
of deformation, as well as shape and environment. Of particular 
importance is the sections on "dynamics" properties of rubber and 
methods of measuring them. 

13 

‘ Gottenberg, W. G. : "Results of Minuteman Engine Associate Contractors 
Round Robin Propellant Mechanical Properties Test", STL Report, 6121-7244- 
KU 000, July 26, 1963. 

Values of stress relaxation and dynamic modulus as measured by 
Aerojet General Corporation, Hercules Powder Co, , Thiokol Chemical 
Corp. , and Space Technology Labs are compared to show the dispersion 
of data. 

14 

’ Kingsbury, H. B. ; Vinson, J. R. j and Soler, A. I. : "On the Lobar 
and Longitudinal Vibrations of Solid Propellant Rocket Motors", AIAA Paper 
No. 65-172. 


Approximate analytical methods are developed to determine the 
lobar (breathing) mode shapes and natural frequencies of a solid 
propellant motor, the axially symmetric propellant mode shapes and 
natural frequencies during longitudinal vibration, and the propellant 
stresses during forced longitudinal vibration. For lobar vibrations, 
general expressions for both kinetic energy and strain energy are 
formulated utilizing a plane strain solution of the motor undergoing 
arbitrary inextensional deformations of the case. Employing a varia- 
tional technique, analytical expressions for the lobar frequencies and 
mode shapes are determined. For longitudinal vibrations, uhe pro- 
pellant is considered as a thick, hollow, elastic, finite cylinder. 
Functional forms for stresses and displacements are assumed with 
arbitrary radial and axial dependence. The repeated use of Reissner's 
Principle yields equations applicable to various boundary conditions. 
From these, natural frequencies and mode shapes are determined. 
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15 . 

'Kruse, R. B. ; "Laboratory Characterization of Solid Propellant 
Mechanical Properties, 11 AIAA Paper, No. 65-147. 

This paper surveys the development of techniques for the laboratory 
testing of solid propellants, reviews the results to date in regard to 
to the degree of understanding they provide about the nature of these 
materials, and indicates briefly the direction of f uture work in this area. 

NASA SP-8064: "Solid Propellant Selection and Characterization" 

The purpose of this monograph is to organize and present, for effective 
use in design, the significant experience and knowledge accumulated 
in development and operational programs to date. It reviews and assesses 
current design practices, and from them establishes firm guidance for 
achieving greater consistency in design, increased reliability in the end 
produce, and greater efficiency in the design effort. 

17 

' Nielsen, L. E. : Mechanical Properties of Polymers , Reinhold 
Publishing Corporation, New York, 1962. 


This text presents a concise review of a wide variety of mechanical 
properties of high polymers from both the theoretical and experimental 
viewpoints. Of particular importance are the discussions on the time 
temperature superposition principle, models of viscoelastic theory, 
and interrelation ships between various properties. 

Stibor, G.S. : Working Papers 


A study was conducted to investigate the structural dynamics of a 
solid rocket motor and attachment structure. The discussion and 
results of the study which follows consists of the following: structural 
resonant frequencies of the motor and motor with attachments; motor 
chamber acoustic modes (self- gene rated vibration environment); and 
vibration data from previous static tests of 156-inch diameter solid 
rocket motors. 


Williams, M. L. : "Structural Analysis of Yiscoelastic Materials, ” 
AIAA J., 2, 735-804, 1964. 


This paper discusses the material characterization, engineering 
analysis, and failure criteria of viscoelastic materials. Of parti- 
cular importance is the material characterization where different 
mathematical models are used to explain the behavior of viscoelastic 
material. 
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20 

Williams, M. L. ; Landel, R. F. ; and Ferry, J. D. ; "The 
Temperature Dependence of Relaxation Mechanisms in Amorphous Polymers 
and Other Glass -Forming Liquids, "J. AM. Chem. Soc. , 77, p. 3701, 1955. 

The temperature dependence of stress relaxation of polymers is 
studied. Data taken at different temperatures can be shifted on 
the time scale to obtain relaxation data over broad time periods. 


Section 2. 

21 ' 

'Achenbach, J. D. : "Dynamic Response of a Long Case-Bonded 

Viscoelastic Cylinder, 11 A3AA Jorunal, Vol. 3, No. 4, pp. 673-677, April 1965. 

22 

"Achenbach, J. D. : "Dynamic Response of an Encased Elastic 
Cylinder with Ablating Inner Surface, 11 AIAA Journal, Vol. 3, No. 6, pp. 1142-1144, 
June 1965. 

22 

Achenbach, J. D. : "Thickness Shear Vibrations of an Ablating 
Rocket, " TR No. 65-3, Northwestern University, July 1965. 

24 

•Achenbach, J. D. : "Forced Vibrations of a Burning Rocket, " 

AIAA Journal , Vol, 3, No. 7, pp. 1333-1336, July 1965. 

•Achenbach, J. D. : "Dynamic Response of a Viscoelastic Cylinder 
with Ablating Inner Surface, " Northwestern University TR No. 65-6, Dec. 1965. 
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for Acoustic Mode Shapes and Frequencies in P«.ocket Motor Combustion Cavities, " 
Unpublished Report, Hercules Incorporated. 

Q A 

’Anon: Aerospace Metals Handbook, Vol. I, "Ferrous Alloys, " 

Code 1213, March 1966. 



Anon, : ''Fracture Toughess (FZM-12-408)" General Dynamics 
Corporation, Fortworth 

32. 

' Baltrukonis, J. H. and Gottenberg, W. G. : "Transverse Wave 
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